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On z-dimensional projectively flat spaces admitting 


a group of affine motions G; of order 7> 7?—n 
By 
Yosio MUTO 


Synopsis 


The necessary and sufficient conditions that an z-dimensional manifold 
A, with symmetric affine connection and with vanishing projective curvature 
admit a group of affine motions G, of order +> —n are obtained in the 
form of tensor equations for z= 3, by studying the integrability conditions of 
the differential equations of affine motions 


Sa Bey aR vd 


We find that we have three types of spaces. The spaces of the first type, 
T1, admit a transitive group G, with r= 7, or an intransitive group G, with 
y=n’—1. The spaces of the second type, T2, admit a transitive group G,, 
7 =n’—n-+1, and the spaces of the third type, T3, also admit a transitive 
group G,, r=n’?—n+1. The rank of the symmetric part of the Ricci tensor is 
one for Tl and T2, while it is two for T3. The skew symmetric part vanishes 
for Tl and T3, while it is of rank two for T2. The existence of such types 
are proved. As, according to G. Vranceanu, a projectively not flat space does 
not allow a group of affine motions of order more than 2—27+5, we find 
that we can omit the words “ with vanishing projective curvature” for m= 5. 


Introduction 


We say that an z-dimensional affinely connected manifold An with 7 as 
parameter of connection allows an infinitesimal affine motion when the Lie 
derivative XJ’, with respect to the vector of infinitesimal transformation 


&* vanishes: 
(1) . AT =0. 


i 


As we consider only symmetric connections, (1) is equivalent to 
(2) A= Raat" 


where 
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(3) cos 


(4) Rie =D le Ups tT way 


* Vo BY, © Bw, v 


the semicolon and the comma indicating covariant and ordinary partial dif- 
ferentiation respectively [K. Yano].” 

The problem whether an A, admits an y-parameter group of affine mo- 
tions G, or not will be solved by studying the integrability conditions of the 
system of partial differential equations (2) and (3). If G, is the complete 
group of affine motions in Ay, the condition of integrability of (2) and (3) is 
obtained as a set of 2?+nu—yr independent homogeneous linear equations in 
é* and", 


(5) EXV7,tiV.,=0 (A=1,--,m; m=n'?+n-7r), 


for the equations (2) and (3) and the equations obtained by differentiating 
these covariantly are all linear and homogeneous in &* and &*. In the pre- 
sent paper a projectively flat A, is studied and the conditions that it admit a 
group Gy of order > n?—un are obtained. 


§1. The four types of projectively flat An 
The curvature affinor of a projectively flat An has the form 
(6) R’ Vw re of Oi RD OD a is 3 


As the integrability condition of (2) and (3) is that the set of equations, 


(7.0) Ah peo, 
C7) EE Sale Bye yf — 
(7. p) CR (0 


where &* and &* are the unknowns, be compatible, we find that this set must 
be equivalent to (5) algebraically [K. Yano]. Hence there can be at most 


1) See references. The Lie derivative of the affine connection Hast is given by 
XThy= eg, yet bot onlay ee Ou Ed +h, v 
and the Lie derivative of an affinor TX is given by 
XT bY Te ga Te a ATR hee 
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o 


only 2 7—1 linearly independent equations in (7.0), for we are studying the 
case 7 > #—n. 


Substituting (6) into (7.0) we get 
(8) A = Qs 


and decomposing //,, into the symmetric part S,, and the skew symmetric 
part V,,. 


(9) Ey gt 
we get 

(10) XS,, =0 
and 

(11) XV, = 0. 


First we consider equations (10). Transposing the terms containing é*to 
the second member we get equations of the form 


(12) Sua E% +Say Es = eae 


As Si is symmetric in » and v, we can find out a coordinate system such 
that at the point in consideration, (%)o , we have 


Suv = €u Suv , & = +1, 0, 
so that (12) becomes 
(13) ep Et +ey EY = Bye c 
If the rank of Su, is more than two, wecan take e, 4-0 for ux=1, 2, 3, and we 
get 3 7—3 independent equations by putting p=1, .,2 and vy =1,2, 3. This 
contradicts the assumption 7 > m?—n, and we find that S,, must satisfy 
(14) Sw =a Sp Sy $67, 7y, Ca,6= 41, 0), 
where S, and Ty are linearly independent. 


11 First, we shall consider the case where a,b = +1. 
Substituting (14) into (10) we get 


(15) a{(XS,) Sy + Sy CX Sy} + 0{CK T, ) Ty +Ty CX Ty $= 0. 


Multiply by s” and constract, where s* satisfies Ss s* =1 and Tas” =0, then we 
shall find that the vector XS, is a linear combination of Si and Tw, 


x Sp = ak Sa Bale 
Similarly we get 
XT Suck Oba. 
Substituting these into (15) we shall find 
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aa=aBt+by=b6=0, 


hence we get 


(16) XSi=B8 Ty, 
(17) XT,= 7S, 
(18) aB+tby=0. 


If @ and y are eliminated, these give just 2”—1 independent equations. 
Since we assume 7 > 22—x there can be no new relation in (11). 

If we choose a coordinate system such that, at the point (x*)o, Su=sl 
and T, = & are satisfied, then the equations (16) and (17) give equations 
expressing £1 and £2 in terms of linear forms of &*, exeept &} and &, which 
obey 2&5 = B+---, £2=y+---, hence 


(19) afl+bi2—-.., 

-- in the second member denoting terms of &«. As there can be no more 
independent equation, the remaining unknowns ¢2 can take arbitrary values 
independently at (x%*)o.2 


If we transpose the terms containing é* to the ‘second member in (11) 
we get 


. (20) Vue EF —Vie EG = 


for V,, is a skew symmetric tensor. If in the first member —V,3, the 
coefficient of £3, does not vanish, £3 can not be independent, contrary to the 
assumption 7 >2?—n. Hence we get V,s =0, and similarly V,,=0. As Vu, is 


skew symmetric, the only non-zero components are V2 and V2. Returning 
to a general coordinate system we get 


(21) Wy ーー V (Sy Ty —S, Tu ) 


which must hold at any point. 

Substituting (21) into (11) we get XV =0, that is, V,4 €*=0 because of 
(16) and (17). This equation is independent of the equations hitherto obtained 
if the gradient V,« does not vanish. Hence we get V = constant. 

1.2 Next, consider the case where b6=0. We have 


(22) | Sw = aS, Sy, a=+1 
and substituting this into (10) we get (XS,)S,+S,(XS,)=0. From this 


relation we soon find 


2) We use the indices in the following manner: A,p,--=1,---, m3 i. jp =2,0°, 3 
a, b, Sachs … 。 2. 
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(23) a5, 0. 


If we choose a coordinate system such that, at a point (x)y, S,= 6! are 
satisfied, et are expressed as linear forms in &*. Then we get 


(24) Vu bi Va £4 = 


from (20), where the second member is a linear form in E*. As we have n 
linearly independent equations in (23), there can be no more than n—1 
independent relations in (24). : 

Now we can find out a coordinate system such that at (x Yo Via = 0 are 
satisfied beside S,= 61. Put w=1 and v=7 in (24). Then if V2<0 we find 
that we can express £2 in terms of e, and &*, 

We have already 2m—1 independent equations, and the remaining equa- 
tions 


Viatg—Via £9 +V jn £2 —Vin EX =---, 


which are obtained from (24) by putting »=j and v=k, must be satisfied 
whenever the former ones are satisfied. In other words, we have two sets 
of equations 


(25) éL=0 (mod 44 &), £&=0 (mod £4 &*) 
and 
(26) Via E4—Vya E4=0 (mod £4, é«), 


and the equations (26) are some linear combinations of the equations of (25). 
(Such equations as (25) and (26) mean that the first members are some 
linear forms of & and &*.) But obviously this can be the case only when 
Vja=0. On the other hand we have already Vig =0, and consequently we 
find that the only non-zero components are Vy and Vy. Returning to a 
general coordinate system we have 


(27) Ves a Su qe, —S, wa. 


where 7, is some covariant vector. 

If Vz, =0, we can not take Vj23<0. In this case we can show that ve 
vanishes. For, if Vuy2<=0 and V;,=0, we can choose a coordinate system such 
that at («*)) we have V23=1, V4 = = Vin = 0. We get from (24) 


Vg f=, 


人 8 =" 
(28) 


Vx £3, —Vni & = 
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by putting =2and v=1, 3,., %. These are evidently m—1 linearly inde- 
pendent equations, and there can be no more relation in (24). But if we put 
p=l1 and v=3,--, 2, we get VaiEi=--, --, Vi i=, which are dependent of 
(28) only when Va; =0, contrary to the hypothesis Vy, +¢ 0. 

Hence we know that we get (27) or Vu, =0 whenever we have (22). 
13 If z=b=0 in (14), we get Rr ,,,=9, fora projectively flat space with 
skew symmetric Ricci tensor is flat [I. P. Egorov, Y. Muto]. 

Collecting the results obtained and renewing some of the letters used we 
get the 

Lemma 1. Jf a curved n-dimensional manifold Ay, with symmetric affine 
connection and with vanishing projective curvature, where n=3, admits a 
group of affine motions Gy with order r > n’—n, it belongs to one of the four 
types of spaces, with curvature affinor satisfying 


T1 RD, = GAy(& Ay — 84 A,), 


* BVO 


T2 R* ayo = GA, (8 A, — 84 A,) 
+28 (A, By —A. By +8 (Ay Bu —Ao By ) 
—& (A, By —A,B,), 

T3 RX 4 = 8 (a Ap Ao +b By B, )—8 (aA; Ay +0 B,B,), 


TANTO) 


T4 R* Bye =o (aA, Aw +6 Bu See )—8 (aA, Ay +5 Bu B, ) 
+ k [2 OX CA, Be —A, BB +e (Au ‘Bs —Aw Bi ) 
ーー (Ay fey —A, B, J ’ 


where a,b= +1, k is a non-zero constant, and A, and B, are linearly inde- 
pendent vectors. 


§2. A, of the first type, T1 
For Tl we get 


(29) XA,=0 


for we had (23). From (29) we get (X A, );,=0, hence X(A, -y )=0 because 
of PD her = 0(K. Yano]. Put A,,=A,; y, then we get 


(30) XK Ay =0, 
that is, 

(30.1) X Aqy =0, 
(30.2) X Aq =0. 


As (30) are equations of the same form as (8), we can deduce that Any 
satisfy one of the following sets of equations: 
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(31) fees 02,1, +mM, My; 
Atuy) =p (Ly MI M.) ’ 
(32) に = LL, Ly, 
Atv} = Lu M, —L, Mi,, 


(33) Ayy = some skew symmetric affinor. 


Here 7 and m are +1 and fp is a constant, while Z, and M, are linearly 
independent vectors in (31). In (32) Z,3<0, but M@, may be dependent of 
L,. Moreover we can get the equations 


(34) XL,=1AM,, XM=pL,, KAKX+Fmp=0 
from (31), and 
(35) Xela 0 


from (32), just as we get (16), (17), and (18) from (14), and (23) from (22). 

If (31) is satisfied and the vectors A,, Ls, and M, are linearly independent, 
we have 3v—1 linearly independent equations in (29) and (34), contrary to 
the assumption 7>n?—n. Hence we can put A,=ql,+7M,. Take a co- 
ordinate system such that at a point (%*)o we have Ly =61 and M, = &, then 
from (34) we get 2n—2 equations expressing £1 and 2 in terms of &* , except 
él and &?, which satisfy 1&}+mé2=0 (mod &*). On the other hand from (29) 
we get gé1 +r £2 =0 (mod &*) which give g£}=7é?=0 (mod &«) if we put 
k=1,;2. We have 2 equations in all except the case g=r=0 which can 
not take place because of Ay*<0. This contradicts r > n’—n. 

Consider the case where (32) are satisfied. If Ls and A, are linearly 
independent, we have 2m independent equations from (29) and (35), contrary 
to r>n’—n. Hence we get L,= XA, and 


(36) Aggy =O Lip Ag 
Applying (32) and (35) to (30.2) we get 
(37) XM = BAz. 
If (33) are satisfied, we get 
(38) A = Ay By —Ay Bu 
from (29) and (30.2), just as we got (27) from (23) and (11). Substituting 
(38) into (30.2) we get 
(39) A Be = 0 Ag 
Renewing some of the letters we can express the results obtained as fol- 


lows: 
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(40) A = Au Ay +An Bi, —Ay, lath a 


If a is zero we get the case of (33). Remember that this is the result ob- 


tained from (29) and r>7?—x. 

Now we can show that A,B,—A,B, in (40) must be zero. Applying 
the expression of the curvature affinor of Tl to the Bianchi’s identity 
Fes =0, we get 


° lve ;k] 
Aj Ay +Ap Avo—Apy Ao —Ap Aw =0, 
from which we get 
3 A, (A, Bu —An By) = 0 
because of (40). Hence we have 
(41) Ay =a ApAy. 


The case a =constant has been completely studied CY. Muto]. If a@ is 
not a constant, we get 


(42) KA=A,pe*=0 

by substituting (41) into (30). Differentiating (42) and applying AT*, =0, we 
get 

(43) ACs Sen 


Now, if the gradient vector a,, is linearly independent of A,, we have 2”+1 
independent equations in (29), (42), and (43), contrary to r>2’—mn. Hence 
we get 


a, kh = B A, ’ 
and substituting this into (43), we get X8=0. This gives 
ioe kh = oY A, ’ 


and from these we get 

(44) Oy =f (a) Ay. 
Then (42) becomes 

(45) Avé*=0, 


and the group G; is given by (29) and (45). 

Hence it is intransitive and 7 = #?—1. 

As a result we have the 

Lemma 2. A space of type Tl may have a group of affine motions G; with 
r > n?—n only when (41) is satisfied. If a is a constant, Gr with r= exists 
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and it is transitive. If a is not a constant it must satisfy (44) and the group, 
if it exists, is determined by (29) and (45), it is intransitive and + = ne—l, 


§ 3. Some preparations 


For an An of type T2 we have (29) and X(A, By —Ay Bz) =0, hence 
(46) XK By = OA 


This equation is the same as (39), but the vector B, is determined by 
Re wie = —(4+1) (A, By —AyB,) , not by (38). From (29) we get (40) again, 
but in order to avoid confusion we must change the letter B for another, for 
example C, hence 


(47) Ap. Au A, +Ay Cs —A, Cr . 


From (47) and (30.2) we get XC, =p Ay, and, if the three vectors A,, Bu 
and C, are linearly independent, we have 3z—2 linearly independent equations, 
contrary to r > n’—n, hence C, must be some linear combination of A, and. 
B., so that we get 


(48) Auy — a Au A, +Y CAS B, —Ay, B.D) . 


We have already 2x—1 equations in (29) and (46), so that the equations 
X a=0 and Xy=0, which are obtained by substituting (29) and (46) into the 
equations obtained by substituting (48) into (30), must be some linear com- 
binations of (29) and (46). But this can not be possible except the case of 
constant a and y, for (29) and (46) are equations expressing &* in terms of 
linear forms of &*, whereas Xa and Xy contain no &*. 

Consequently we know that a necessary condition for r>mn?—n is (48). 
with constant a« and y, and the group must satisfy 


(49) MAg = 0,me = BAy.. 
For an A, of type T3 or T4 we have a,b+-0 in (14). We get 
(50) XA,=cB Bu, 
(51) XB, =B8 A,, 
with 
(52) —c=b/a= +1, 


just as we got (16), (17) and (18) from (10) and (14). 
We can consider that (49) is a special case of (50) and (51) if we allow 
c =0, and we shall treat these equations for a while. 
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Take such a coordinate system that at the point in consideration we have 
A, = él and, Be= oe Then (50) and (51) become 
El + Aya E*=c BO £24 Bus E* =f 6h 


fe? 


at this point, so that the 2n—1 quantities 


are obtained as linear forms of &*®. Because of x >m?—n no other relation 
is possible, and £2 and &? are arbitrary. 
Differentiating (51) covariantly we get 


(54) X Buy—PBy Au—B Ap, = 0 


because of XJ". =0. Substituting 8 =&? (mod &* ), which is obtained by put- 
ting “#=1 in (51), into (54) and putting w=7, we get 


(55) Bja E94 + Ba, E% —Ajyé?=0 (mod & ). 
Putting j=a and v=6 in these equations we get 
Bac &¢ + Bes ta —Aap &2 =0 (mod &*), 


for the quantities in (53) can be crossed out. As ¢&c and £? are arbitrary, 
their coefficients must be zero, hence 


(56) Ag=0, Ba = 
If we put j7=a and v = 1 in (55), we get 
Baz 3+ Be £5 —Aqi £2 = 0 (mod &* ) 
because of (56), and we know that 
(57) Ba =O)" Bore An: 
If we put 7=a and v =2 in (55), we get 
Bez 5 —Aaz £2 = 0 (mod &* ) 


because of (56) and (57), and we know that Baz=Aa.=0. These results are 
summarized in 


(58) Aay = Ba,=0. 


3) In the following Any, Buy, Apvo ee inea Au: v b Bu: y, Ap; Veo. re 
4) See foot-note (2) in page 4. 
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Moreover, if we put 7=2 andv=4, 


we get 

Bye £5 + Bys £1—Avs £2 = 0 (mod £*), 
hence 
(59) Bo =0, Ax =c By, 


and, if we put 7=2 and vy = 1, we get 

Boz 2+ By EL—Ay £2? =0 (mod é*), 
hence 
(60) Ay = Bytc By. 
Finally, if we put 7=2 and v = 2, we get 

By + By ELA £2? =0 (mod E*), 
hence 
(61) Ax = ¢ (By2+ By) . 


A similar deduction is possible about (50), but the same result is obtained 
by exchanging A and B, 1 and 2, ¢ and 1/c in the equations obtained above, 
so that the whole result is 

Aa = Bay = Aja = Ba =0, 
A = c¢ Bia, 
(62) Ay =c By+ Bz, 
Az = c (By2+ Ba), 
c Ay + Ax =c By, 
Ayp+Ag, =c By. 


Now we can return to a general coordinate system, and get 


Apy = Gy, Ay Ay + dhe (A, B, —A, Bu), 
(63) 
Buy = Au Cy +52, Bu Ay —ay2 By By 


when ¢ is zero, and 


Auy — 1 Ay Ay +m (A, By —Ay Bude LB, By +c BuCy ’ 


(64) 
Bus = =/ Bu Ay —m BB, +A, GS 
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when c= +1, for we get Ajy = —B and As; = —By from (62). C, is some 
covariant vector. These equations are to be satisfied throughout the space. 


§4. A, of type T2 


As we have already (48) for T2, we get 
(65) Ayy = a Ay Ay +y (Ay By —A, Bu), 
(66) Buy = Ay CG; +6 B, Ay ay Bu ep 


from (63), where a and y are constant. Differentiating (65) covariantly and 
substituting (65) and (66) into the resulting equations, we get after some 
calculation 


Auyvo—Ayoy = (2a—8) y A, (A, Bo —Aaw By). 
‘On the other hand we have 
Ayyo—Apoy = —3 Ay (Ay Bu—Aw B, ) 
because of the Ricci identity, Apyo—Ayov = —R ya 
(67) 2ay—sy = —3 


so that 6 is a constant. 
From (66) we can derive 


Aa. Hence we get 


Buyo— Buoy — Ay (Cyw—Cay )—Ca—8) Ay (A, Ge —A.w Cs ) 
—27yA,(B, Co—B, C, )+38yB, (Ay Bs-AB,) 


where Cyo = Cy: o « But we have 
Byvo—Buoy = —R* ve Ba 
= (aA,—3 B,) (Ay Bu—Aw By), 


hence we get 


(68) Cy Cw—(e—8) (A, Co—A, C,) 
ー2 ry, (By, Cu —Bu Cy)—a (Ay ibs —A.w By — (0) 4 
(69) $y=—I1. 


Now, for T2 we have 
(70) igs =a A, Ay, トス A B, —Ay, Bs . 
‘On the other hand we know that we can derive the important relation 


(71) Huy yo TF zn sy = 0 
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from (6) for x > 2. Substituting (70) into (71) we get 
(aA, Ay +A, B, Ay B,): o 

—(@.A, As+A, Bo Ap Bp! y= 0: 

These equations are reduced to 
3 ay—a—éd = 0 

because of (65) and (66), so that we get 
(72) ayye=—1, a=—2/y, 6=—I1/y7, 


because of (67) and (69). If we consider only real numbers we get a = —1. 
If we differentiate (68) covariantly and take the skew symmetric part of 
the resulting expression, we get no new relation, for we have Crwai=0 and 
a—26=0. 
After all we get the 


LEMMA 3. A space of type T2 may have a group of affine motions Gy with 
r> n—n only when A, and B, satisfy (65) and (66), where a, y, and 8 are 
determined by (72) and C, satisfy (68). We get r=n’—n+1 and the group, 
if it exists, is determined by (49), hence transitive. Besides, if we consider 
only veal numbers, a must be —1. 

An example of such an A, is introduced by I. P. Egorov [I. P. Egorov]. 


§5. A, of type T3 and type T4 
We get (64) for T3, where c= —b/a. Substituting 
fa = aA, A, $08, B, 
into (71) we get 
a do Ay FA, A A Ao A A 
+6 (Bio By +B, Bs—Bis Bo—B, Bov) =0, 
and substituting (64) into this we find 
{3amA,+(2bl—acl) B,} (A, Bu —AuB, ) 
+(ac+b) B, (Ay C.—AwC, )+(ac+b) A, (By Co—BoC, )=0, 


so that =m=0, for A, and B, are linearly independent and act+b=0. Hence 
(64) becomes 


(73) Ai n=6 By Cy 5 Buy = A, G; . 
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For T4 we get 
a(Ass Ay +A, Avo—Auw Av —Ajp Aov) 
+6 (Buo By +B, Byo—Buy Bo —B, Boy) 
+R (Avo By +Ay Bio— Avo By— Ay Buw 
—Ay B.—A, Boyt Avy By + AaB) = 0 
from (71) and 
Tuy = a Ay Ay +6 B, B, +k (Ay By —A, By). 


As the third term in the first member vanishes because of (64), we find (73) 
for T4, too. 
Differentiating (73) covariantly we get 


Aine = CA Gy Ce +c Bu Cus ’ 
(74) 
Base =C Bu Cc, Gs Au Cs , 


and substituting these equations into the Ricci identities we find 
¢ Bu (Cyw—Cayv) = —b By (Ay Bu —AwB, ), 
Ay (Cywo—Cov) = @ Au (Ay Bo —Aw By ) 
for T3 and 
¢ Bu (Cyw—Cov) = —b By (Ay Bo —Aw B, )—3 RAy (Ay Bo —Aw By ), 
Ay (Cyo—Coy) = a A, (A, B, —Ao By )—3 k Bs (Ay Bu —Aw By ) 
for T4, hence 
(75) Cyw—Cay = &@ (Ay Bu —An By ) 


and k= 0. 

Consequently we get the 

LemMa 4. A space of type T3 may have a group of affine motions G, with 
r>n—n only when A, and B, satisfy (73), where c= —b/a and C, satisfy 
(75). The group, if it exists, is determined by (50) and (51), hence transitive, 
and r=n’—n+1. For spaces of type T4 7 can not be more than n2—n. 


§6. Necessary and sufficient conditions 


Now we prove the 
THEOREM: The necessary and sufficient condition that a curved projec- 
tively flat n-dimensional manifold Ayn with symmetric affine connection admit 
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a complete group of affine motions Gy of order r > n2—n is that the curvature 
affinor belong to one of the following three types T1, T2 and T3, where the 
vectors satisfy the adjunct equations. Such connections and the groups really 
exist. 

The curvature affinor 


(/6.1)— 11: k* = a A, (oO Ap —O* Ay ), 
C= Tore 2205 
Rid) TZ RE a Ay (0% A, 0 A; ) 
+20 (A, B, —Aw By )+& (Ay Be —Aa B, ) 
SUSCAER RD BE : 
a=+1, A, and B, are linearly independent, 


03.4) 913 Re 


Sve = 2 Ay (0% Ag—& Ay TO B, (8 By —S B, ), 
a,b=+1, A, and B, are linearly independent. 


The adjunct equations 


for TL 

(76.2) Ag; va aA. A, A 

(76.3) a, v=f (a) Ay, 

FOr 

WAD C0 A, A, $A B, —A, B,D) Y 
Buly=AuCy+ay Ay, By —7 Bu By, 

(77.3) (En Hag) —C, Faye — ay (Az C, —Ay, Cm) 

+2y (B,C, —By, C, )+a CA, B, —A, Bu), 
ay?=-—1, 

for T3 

(78.2) fi ;v =—ab BCs, 
Busy = AuCr, 

(78.3) Cp 3 —C, Pye =a (Au By —A, Bu ) . 

The groups 

for LL 

(76.4) XAu=0; Ka =08 


r= 72, transitive, if @ is constant, 
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7 = n?—1, intransitive, if @ is not constant, 


for T2 

(77.4) X ADO; 0B = Ap’, 
7¥=nv—n-+l1, transitive, 

for T3 

(78.4) X A, = —abBP Bi AB, =B Az, 


r=n’—n+I1, transitive. 


In these equations (77.2) and (77.3) are direct results of (65), (66), (68) 
and (72). £ in (77.4) and (78.4) may be eliminated or we may add equations 


B,y=XC,—ayB A,—278 B, 
for T2, and 
B, y= 7a, 


for T3, and take 8 as an unknown function beside &* and —*. These equa- 
tions are obtained by differentiating XB, = @ A, covariantly and taking into 
account of CXB,) :y= X (By: v), (772), and (78.2). 
We need only to prove the existence of the connections and the groups, 
for we have shown that the condition is necessary. The proof is given short. 
For Tl the existence of connection is proved if the following mixed 


system of partial differential equations for the unknowns /”%,, A,, and a 
(79.1) Liat, 
(79.2) Ley ot pate ek ie keto dot GA, (ONAN A 
(79.3) ee v= Tt, Aata Ay Ay, 

ay = Ff (a) Ay 


admits solutions. But the equations obtained by differentiating (79.1) partial- 
ly give R* -.,,.;=0, which is satisfied because of (76.1), that is, (79.2). The 
equations obtained by differentiating (79.2) partially give RX tee = 9, 
which is satisfied because of (76.1) and (76.2), that is, (79.3). Finally the 
equations obtained by differentiating (79.3) partially give relations which are 
satisfied because of (79.2) and (79.3). Hence (79) admits solutions and the 
connection exists. 
The existence of a group is proved if the mixed system 

(80.1) B= 


My 


(80.2) ey Py sit Eve 


EE a  EEEEEEEoEEEEEeEEeEeEeEeEeEee—eEeeeeee eee 
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(80.3) XAn=0, Masg 


admits solutions é* and &*. But the equations obtained by differentiating 
(80.1) covariantly give nothing more than (80.2), for we have CO meer Cees 
See = RR Bon ae wae’. The equations obtained by differentiating 
(80.2) covariantly give nothing more than X R* we=0, Which is fulfilled be- 
cause of (80.3). Differentiating (80.3) covariantly we get X (Ay; 1) =0 and 
X(a.,)=0, for (80.2) is equivalent to XJ ee = 0. But these are satisfied 
because of (79.3) and (80.3). Hence the system (80) is completely integrable,. 
and we can get a transitive group G,, r= 2, if a is constant and an intransi- 
tive group G,, r= n?—1, if @ is not constant. 

There is no special difficulty in proving the existence of connections for 
4 WAR tes ta 

For T2 we have 


(81.1) Led BeOS 
(81.2) BS ent AG eo YL 


+a A, (8 A, —8% A, )+28 (Ay Bua—A. By ) 
+8 (A, Be —Aw B, )—8& (A, B, —Ay By), 


A, yal %, Aat(2a A, A, +A, By —Ay Bu) y, 


(81.3) 
oe = Ue Bs+A,C, +a y Bu Ay may BwBy, 


(81.4) Cy, v—Cy, uw —ay (Au Gs —Ay Gi +2 FY, (Bu C, —B, Cu ) 
+a (A, B, —Ay Bud . 


But the equations obtained by differentiating (81.1) partially give equations 
fulfilled by (81.2), the equations obtained by differentiating (81.2) and (81.3) 
partially give equations fulfilled by (81.2), (81.3) and (81.4), and finally the 
eguations obtained by differentiating (81.4) partially give equations fulfilled 
by (81.3) and (81.4), as we have ay2=—1. These facts are already shown in 
§ 4 and the connection exists. 

For T3 proof is got in like manner. 

To prove the existence of a group for T2 we consider the mixed system. 


(82.1) = om 
(82.2) i sv = ere ve ee ’ 
(82.3) Ave 0 XB AA 


(82.4) B,v=XC,—ayBA,—2y8B,, 
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where &*, £*, and 8 are the unknowns. The equations obtained by differen- 
tiating (82.3) covariantly are fulfilled because of (81.3), (82:3) and (82.4), 
while the curl derived from the second member of (82.4) vanishes because of 
(81.3), (81.4), (82.3) and (82.4). Hence (82) is completely integrable and we 
can get a transitive group G,, r= n?—n+1. 

For T3 proof is got in like manner. 


In concluding the paper, the author wishes to express his hearty thanks 
to Professor K. Yano for his kind encouragement and suggestions. 
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Viscoelastic Properties of Polymeric Substances* 
By 


Keiichi SHIBATA 


Abstract 


(1) Measurements on the transient behaviour of polymethyl methacrylate 
reveal that the torsional and bending creeps are linear with logarithmic time 
over a considerable time range. 

(2) Comparison between the creep data and theory of Tobolsky and 
Eyring shows that somewhat ambiguous modifications have to be made for 
the parameters in their fundamental equation. 

(3) By means of Gross’ theory, a distribution function of retardation 
times is deduced from the creep data. 

(4) Measurements of dynamic Young’s modulus of polymethyl metha- 
crylate are made with bending vibration frequencies of 500 and 100 cycles 
and at various temperature ranging from a room temperature up to 120°C. 
A Young’s modulus-temperature diagram constructed from the data of several 
investigators including the present ones shows that there is an anomalous 
region in the vicinity of 100°C where the Young’s modulus is markedly 
dependent upon temperature and frequency. An attempt is made to correlate 
distribution function of retardation times obtained from the creep data and 
the vibrational experiment. ) 

(5) Correlations between the transient and the steady-state behaviours 
suggested by Dunell and Tobolsky are discussed on the basis of Gross’ theory. 
If the stress relaxation (or creep) is found to be linear with logarithmic time 
over a considerable time range and the principle of superposition is valid, the 
simple correlation between the rate of stress relaxation (or creep) and the 
elastic loss Cor imaginary component of compliance) is derived so long as 
the relaxation Cor retardation) frequency function is of the rectangular form 
and the distribution function of relaxation (or retardation) times is of the 
form of one type of Pearson’s distribution function. The experimental veri- 


fication is made. 


* The researches were carried out at Kinoshita Laboratory, Scientific Research Insti- 
tute, Tokyo, from 1948 to 1952. 
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I. Introduction 


In recent years, the great development in the production technique of 
high polymers has brought a renewed interest in the studies of their me- 
chanical behaviours. These substances have properties which place them be- 
tween ideal solids and ideal viscous liquids, i.e., they have rigidity as well as 
viscosity. Some attempts have been made to develop a general theory which 
would describe the mechanical behaviour of bodies which exhibit simultaneous 
elastic and viscous behaviours. Although these theories have thrown much 
useful light on many aspects of general problem, various essential behaviours 
‘still remain so far unanalyzed on account of the scarcity of available experi- 
mental data. 

In an experimental study on the mechanical behaviour of viscoelastic 
‘body there are two ways of attack; they are: (a) The measurement of the 
transient phenomena caused by sudden application of a constant load (creep) 
or deformation (stress relaxation). (b) The measurement of the steady- 
state behaviour under alternating load or deformation. Observations are 
made on the transient behaviours for the duration of about one minute or 
longer; the steady-state behaviour involves measurements with frequencies 
higher than 1c/s. 

The present investigation was undertaken in order to carry out 

(a) measuements of torsional as well as bending creep and dynamic 
‘Young’s modulus of polymethyl methacrylate, 

(b) comparison between the experimental results and theories, 

(c) deduction of a distribution function of retardation times from the 
experimental results 
and 

(d) derivation of a formula which correlates the transient and the 
steady-state behaviours. 


II. Experimental (Transient Behaviour) 


Previous experimental researches? have laid the groundwork for an 
investigation of the transient behaviours of polymeric substances. One of the 
main purposes of the present measurement is to acquire information that 
may conduce to exact features of the time dependence of recoverable retarded 
‘deformation (so-called “reversible relaxational strain”). 

By use of torsional deformation of a cylindrical rod as strain, it may be 


1) For instance; Kobeko, Kuvshinsky, Gurevich: Bull. Acad. Sci. U.R.S. 5.5 ou 1937), 
336; S.D. Gehman: J. Appl. Phys., 19 (1948), 456; A. Tobolsky, H. Eyring: J. Chem. Phys., 
11 (1948), 125; H.A. Robinson, R. Ruggy, E. Slantz: J. Appl. Phys., 15 (1944), 348; J. A. 
Sauer, J. Marin, C.C. Hsiao: J, Appl. Phys., 20 (1949), 507. 
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possible to make the measurement with high accuracy, shape and dimension 
of the test piece being unchanged during the creep, and to obtain facilities 
in making the direct comparison between the experimental results and the 
phenomenological theories. As the material to be used in the present investi- 
gation, polymethyl methacrylate™ was selected, because its molecular structure 
is simple and its experimental data are lacking. 

General arrangement of the 
apparatus is shown in Fig. 1, in 
which J is the cylindrical test 
piece, its diameter being 8.50 mm. 
Two shafts arranged in strict 
alignment are firmly coupled to 
it through chucks 2 and 3, and 
are clamped toa stand with vices 
4 and 5, one of the shafts resting 
through a ball bearing 6. A wheel 
having a sufficiently small moment Fig. 1. 
of inertia, its diameter being 
13.6 cm, is attached to the shaft. By hanging a weight 8 over the wheel, the 
test piece is twisted. With the aid of a small mirror 9 affixed at a desired 
position on the test piece, the amount of strain is measured by the method 
of telescope and scale. 

After ascertaining these was no observable strain due to setting of the 
test piece by leaving it unloaded for a week, a load was applied and the 
measurements were made at a room temperature. 

A curve typical of those obtained experi- 
mentally is shown schematically in Fig. 2. On 


sudden application of a load, there is an instan- 3 

taneous elastic deformation, the amount of which a 

is approximately proportional to that of the 7 

load. This is followed by a retarded defor- 3 Time 
mation (creep) which was, in the present experi- Fig, 2. 


ment, still observable after 104 minutes. On 

sudden removal of the load at any desired time during the creep, an instan- 
taneous partial recovery takes place, followed by a retarded recovery (creep 
recovery). There was no observable irreversible deformation in any of the 
present experiments. The instantaneous deformations due to the application 
and the removal of the load seem to be of equal amount,® and the rigidity 
calculated from them gives the values of 1.2~1.1x10!9 dyne/cm?. Attention 


2) Manufactured by the Fuji Kasei & Co., Ltd., Tokyo. 
3) The latter seems to be slightly smaller than the former. 
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should be paid to the fact, however, that it takes a few seconds at least to 
set out the measuement after the application or the removal of a load 
whereby the instantaneous deformation mentioned above (apparent instan- 
taneous deformation) should involve a portion of the retarded deformations 
corresponding to the time interval of a few seconds. The exact value of the 
instantaneous rigidity should be computed from experimental results obtained 
in supersonic frequency range. (The problem is not treated here.) 

Some of the observed values are plotted as shown in Figs. 3, 4 and 5 


a 


Oc =1.6, X107C.8.8.(20°C) 2) end tS i 
ーー - 


ララ 


— ¢ 


/ /0 700 700 
Time (min) ; {nee 


Fig. 3. Creep recovery curve (1/) corresponds to creep curve (1). 
(Unloaded after 2280 min-loading.) 

Creep recovery curve (2/) corresponds to creep curve (2). 
(Unloaded after 10080 min-loading.) 


taking logarithm of time as abscissas. 
From the figures, it may be noted that 
the following equations are quite well 
fitted for the description of the be- 
haviours over a considerable time range 
(1~10! min): 

In the case of the application of a 
load 

y=a+b oy Int, (1) 


and in the case of the removal of the 
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Joad a 
=a!— : °C) 
y=a'—b ap Int (2) quae si2s 


-3 0-3 ) 
where od oo 


y: strain observed, 
oo: constant stress, 
t: time, 
a,a', b: constants independent 
of time. 

It may be expected that 
Eqs. (1) and (2) would not 
describe the behaviour for 
extremely small as well as 
large time scale. 

Measurements were made 
at a room temperature on the 
time dependence of the de- 1 
pression of middle point of 
a rectangular beam under a 
constant load applied at that =a (a) 
point, both ends of the beam 0 : 


70m 
Time (min) 


Strain 


being supported on knife edges. 
Recoverable creep in the Fig. 5. Creep recovery curves (1/), (2/) and (3) 
bending deformation y of the correspond respectively to creep curves (1), (2) 
beam is also described by the and (3). (Unloaded after 10 min-loading.) 
equation 

y=a+t5b' Int (1/) 
which is similar to Eq. (1) over a considerable time range from a few seconds 
up to 3x103 min. The apparent instantaneous Young’s moduli computed from 
the corresponding deformations provide the values of 3.3x10!0 dyne/cm? 
(16°C) and 2.910 dyne/cm? (28°C). 


III. Comparison between Theories and the Experimental 
Results on Creep 


The theories which account for the linear behaviour of creep with 
logarithmic time under constant load are classified into two groups; they are: 
(a) Theory of Tobolsky and Eyring in which Newtonian viscosity in the 
dashpot is replaced by Non-Newtonian viscosity (strictly speaking, Eyring’s 
viscosity), and therefore the solution of the equation does not satisfy the 
principle of superposition. (b) Description of the behaviour by use of an 


4) A. Tobolsky, H. Eyring: J. Chem. Phys., 11 (1943), 125. 
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infinite number of Voigt elements placed in series which is characterized by 
a distribution function of retardation times deduced from experimental data. 
The procedure is entirely based on the assumption that the principle of 
superposition is valid. 


Part 1. Comparison between Theory of Tobolsky and Eyring and the 
Experimental Results. 


Tobolsky and Eyring presented a molecular model to interpret the elastic 
viscous properties of polymeric substances, i.e., they described creep phe- 
nomena in terms of molecular rate and flow processes by using the postu- 
lation that the total rate of strain for the z-th type of unit process is given by 
the following equation : 


dy os Vi: fate BL. nd EN 
We = Gr dhe + tM — exp (FF ) 2 sinh oh ar (3) 


oi: total stress distributed over the unit process, 
mj: number of the unit processes in series per unit length along a 
line perpendicular to a plane over which the stress distributes 
itself, 
N;: number of the unit processes per unit area on the plane over 
which the stress distributes itself, 
G;: appropriate partial modulus to be associated with the unit 
process, 
Aj: average distance projected in the direction of stress between 
equilibrium positions in the relaxation process, 
4F;*: free energy of activiation for the process, 
k: Boltzmann’s constant, 
h: Planck’s constant, 
T: absolute temperature. 


For small loads, it may be considered that the total stress is distributed 
over the primary bond in the molecular model (purely elastic) as well as the 
secondary bond (viscoelastic). In this case the appropriate equations are of 
the forms: 


dy _ 1 doy | 
hore abi? a 
d Y Ye : 

“fe i = + Az sinh Bz a2, . (5) 


5) B. Gross: J. Appl. Phys. 18 (1947), 212; 19 (1948), 257. 
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where 


i) 


kT 
Ap=2 n2r2 pj eXP (= RT (6) 


be gee 
2 2NkT ° (7) 


In the case of application of a constant load, Eqs. (4) and (5) are solved 
under the condition that o>=o1+0.=constant, giving 


nD ae) he Yo) pat (8) 
where 
_ A2BrGi G2 
aes G,+ G2 (9) 


Yo: the value of y at f=0. 
Within the range of 


1 > e7-Bx(00-G SY e- Baloo- Gio) (10) 


one gets an approximate form: 


Se ee ae DS 


On removal of the limitation of 


€- Bz (To-GiY) > e—B2 (@o-Gi Yo) 


y is expressed by the following equation: 


CG BG 2 BG POO (12) 


where 


9 = e~B2 (oo-Gi%0) 
a 

It may be expected from Eq. (12) that the straight line representing the 
relation between y and Inf at medium time scale would bend upward against 
time axis in short time region. On account of the absence of this tendency 
in the present experiment, it seems reasonable to assume that the value of ? 
would be much smaller than a few seconds. 

Supposing that the experimental formula (1) corresponds to the theoretical 


Eq. (11), one gets 


6) Eq. (10) involves the quantities the value of which can not be determined by 
experiment. The range of validity of Eq. (11) is somewhat ambiguous, 
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0 AL 1 2 a 13 
Oi Gi Bo Gi In ve 8 ( ) 
0 BoG, < C ) 


In order to obtain an equation which is valid for the creep recovery, Eqs. 
<4) and (5) are solved under the condition that o;+0,=0. The solution is of 
the form: 


i BiG ppp GD a (15) 


where 
y(0): initial valué of the creep recovery. 
Within the range of 


a >» e~B2Giy >» e7 Bs Gi ¥ (0) 


an approximate expression of 


1 


Bei a 1 (7) 
v= B, com Int (16) 


2 BG 
is obtained. Assuming that Eq. (2) corresponds to Eq. (16), one gets 


if a 


a'=— In Oye (17) 
ee (18) 
BG 


In the present experiment, the load was removed in the stage where the 
creep was still taking place. In view of the experimental results that the 
slopes of the creep and the corresponding creep recovery curves plotted 
against logarithmic time are of the same absolute value, it may be expected 
that the straight lines in y-lné diagram obtained on the removal of load at 
any desired stage of the creep under the same load should be parallel to 
each other. The speculation is verified by the observed values shown in 
Fig. 6. 

The time 4 when y is equal to zero in Eq. (2) is given by Intp=a'/bop. 
By substituting Eqs. (17) and (18) into this equation, one gets 4=2/a. 

Remembering that the value of t depends on the time of removal of 
load, other conditions being the same (cf. Figs. 5 and 6), one may conclude 
that a@ in the case of creep recovery depends upon the time of unloading, 
while @ in the case of creep at a constant temperature is characterized by 
the specific properties of the material. 


7) cf. foot-note (6). 
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Strain 


Time (min) 


Fig. 6. Creep recovery curves (1) and (2) correspond to the same creep curve. 
(1): Unloaded after 90 min-loading. 
(2): Unloaded after 10 min-loading. 


The results of the theory of Tobolsky and Eyring, i.e., Eqs. (11) and 
(16) show that the absolute values of the slopes of creep and the corre- 
sponding creep recovery curves plotted against logarithmic time are not only 
of the same value, but are determined by the characteristic constants of the 
material. The experimental results support the former, while they deny the 
latter showing that the slopes are proportional to the amount of a». Tobolsky 
and Eyring have already pointed out that B, is inversely proportional to the 
amount of go. Hence, it seems reasonable to corclude that the duality in the 
physical meaning of @ mentioned above, originates in that of A». It has 
been left unanalyzed, however, which of the parameters in Az will be re- 
sponsible for the discrepancy. 

In short, in order to interpret. the experimental results by the theory of 
Tobolsky and Eyring, somewhat ambiguous modifications have to be made 
for the parameters A, and B, involved in the fundamental equations. 
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Part 2. Deduction of a Distribution Function of Retardation Times from 
the Experimental Results. 


Some attempts have been made to find out a set of appropriate retar- 
dation times in terms of which experimental results on creep would be described. 
For example, Robinson, Ruggy and Slantz have found three retardation 
times of polymethyl methacrylate with the aid of graphical analysis of their 
tensile creep curves. However, on account of the uncertainity in drawing 
tangents to the creep curves, their graphical method gives its ground to the 
more elegant one presented recently by Gross, i.e., a mathematical deduction 
of an appropriate distribution function of retardation times from experimental 
creep data by considering them as a continuous spectrum. For the determi- 
nation of the complete form of the distribution function, measurement should 
be made over extremely broad time range extending to short as well as to 
long time scale, but it is hard to realize. However, an approximate distri- 
bution function deduced from a set of creep data for considerable time range 
may still serve to discuss qualitatively the mechanical behaviours of high 
polymers. Furthermore, when strain in creep is expressed by a simple 
function of time, the general method of Gross can be replaced by a simpler 
one described below, so long as the distribution function is concerned. 

Assuming that the principle of superposition is valid, the behaviour of 
recoverable creep is described by use of an infinite number of Voigt elements 
arranged in series, i.e., the creep function W(t) is given by a normalized 
distribution function F(T) of retardation times 7 as follows, in terms of 
which the continuous system of Voigt elements is characterized. 


v ()=6|1—| FC eth dr | (19) 


where @ is a normalizing constant. 
The rate of creep p(t) is given as follows by using the retardation 


frequency function p(s): 
p= OF = jes ) e-sds (20) 
where 


| s= 1/;: retardation frequency, 


le F (1) dt =—p(s) ‘ 
As the deformation was recoverable in the present measurement, it seems 


8) H.A. Robinson, R. Ruggy, E. Slantz: J. Appl. Phys., 15 (1944), 343. 
9) B. Gross: J. Appl. Phys., 18 (1947), 212; 19 (1948), 257. 
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reasonable to assume that the amount of strain caused by creep should 
approach asymptotically to a finite equilibrium value with time. To satisfy 
this condition as well as to be able to represent the behaviour of strain in 
short time region, the experimental formula (1) may be modified by intro- 
ducing two parameters Y and ¢ as follows: 


y@)=a"'+ B' In (0 +4)—B' In (€4+2) CAD) 
where 


ae ‘ 
b'=boo , 


and @ and € are constants with the dimension of time, and are much smaller 
than a few seconds and greater than 10! min respectively at a room 
temperature. 

Eq. (21) can be rewritten in the following form: 


ot Svea 29 
y= rl1+ a (22) 
where 
yo=a”+ by in 


Hence, the rate of creep Y(t) in the experimental formula (22) is ex- 
pressed by 


SER OP. SES) 23 
I es aa coe 


By substituting Eq. (23) into Eq. (20) and applying Laplace inverse 
transformation, one gets 


+10 b/ J i 1 ae Ap gas Oe $S_= 9- gs), 
p= 272 eg Yo ie ae Yo 


Therefore, 


/ 9/r —p-S/t 
ARCs ; {° g Nar 


0 T 


where 


eo Tae a Git / 
pial | EO TE 0) 
0 Yo 
Hence, 


Fr)dr= ee) dr (24) 


TN 
In (&/Y) T 
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In the case of ordinary creep where the extremely large value of 7 is 
not taken into account, an approximate form shown kelow may be used, in 
which the term e-$” in Eq. (24) is neglected, i.e., 


F (rt) dr=constant x * e-9't dr (25) 


Eq. (25) is of an unsymmetrical form, its mode being at7=v. (One type 
of Pearson’s distribution function). 

The distribution function (25) has also been given by Voglis as an 
experimental result of dielectric after effect on mica and wax. It will be 
worth making similar measurements on polymethyl methacrylate. 


IV. Experimental (Steady-State Behaviour) 


Measurements of dynamic Young’s modulus of polymethyl methacrylate 
were made with bending vibration frequencies of 500 and 100 cycles and at 
various temperatures ranging from a room temrerature up to 120°C. In 
accordance with the temperature range where the experiment was carried out, 
two methods were employed: one was an electromagnetic method which 
could be applied up to 90°C, the other was an electrostatic one used in the 
temperature region above 90°C. 

For excitation of the bending vibration, an electromagnetic device was 
used, which consisted of a small iron plate cemented on the surface of the 
test piece and a magnetic telephone receiver facing to it at an appropriate 
distance. The oscillating current was picked-up with another electromagnetic 
device of the same construction at the free end of the test piece. It was 
amplified and rectified, and was read on microammeter scale. The mass of 
iron plate was made negligibly small compared with that of the test piece. 
In the present measurements, error due to the mass of iron plate (about 
0.05 g) is estimated to be 1~2%. Correction for rotational energy is neglected 
without significant error. 

In the experiments made in the temperataure region above 90°C, the 
bending vibrations of the beam were excited and picked-up electrostatically 
by using copper foils and plate electrodes in the places of iron plates and 
telephone receivers respectively. Measurement of minute capacity chang of 
the receiving condenser due to the vibration of the beam was made with 
Foster-Seeley’s circuit. 

With either experimental arrangement, a resonance curve was obtained 
by reading the out-put current at various oscillator frequencies. Dynamic 
Young’s modulus £ and the longitudinal viscosity are calculated from the 


10) G.M. Voglis: Zeits. f. Phys. 109 (1938), 52. 
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resonance frequency f, where the out-put current exhibits a maximum value, 
ani from the half value breadth of the resonance curve respectively. (The 
viscosity problem is not treated here.) 

The relation between E and f is shown below: 

B=4m f? (26) 
where 
Z: length of the beam, 
p: density of the beam, 
«x: radius of gyration of the rectangular cross section of the beam, 
kn: constant for each mode of vibrations. 

The fundamental vibration as well as the first and the second modes of 
vibrations were used in the present experiments, the measurements of their 
frequencies being made with Hay’s bridge or Campbell’s bridge. 

In carrying out the measurement, the change of the resonance frequencies 
of the beam of a constant length caused by the elevation of temperature was 
traced. On account of the fact that Young’s modulus, accordingly the reso- 
nance frequency, decreases with increasing temperature, the tracing of 
temperature dependence of Young’s modulus at a definite frequency is not 
easy. Hence, the determination of the value of Young’s modulus at 500 
cycles, for example, was made by interpolating two values obtained at the 
frequencies of the first and the second modes of vibrations, an appropriate 
length of the beam teing selected so that the former was lower than and 
the latter higher than 500 cycles. 

The observed values obtained in the present experiments are plotted as 
shown in Fig. 7 together with those given in literatures. Curves I and II are 
the results of the present measurements in which the former was obtained 
at 500 cycles and the latter at about 100 cycles. The observed values at 2C0 
cycles seem almost fall in with those at 500 cycles in the temperature region 
lower than 90°C. The points 1 and 2 show the values of the apparent 
instantaneous Young’s moduli. (cf. p. 23) Curves A, B and C represent the 
temperature dependence of static Young’s moduli, for which little of the 
experimental conditions was reported. Curve A is the results of measurement 
by Bartoe™ using “quick acting testing machine”. Curves B and C show 
the results obtained by Robinson, Ruggy and Slantz,12 with a thin rod as 
the test piece. Curve B shows the value calculated from the total elongation 
under constant load, and curve C shows that from the apparent instantaneous 
elongation. Hence, it may be presumed that the points I and 2 would lie 


11) Cited in the paper of Robinson, Ruggy and Slantz. (loc. cit.) 
12) loc. cit. 
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on the extension of curve C. The points D and K are the experimental data 
given by Weber and Goeder, the former being calculated from the apparent 
instantaneous elongation as well as compression amount of a rod due to 
loading, and the latter being computed from the compressional vibration 
excited by Kundt’s method or a similar method mentioned in this section. 
It is seen that the data obtained by Weber and Goeder show good agreement 
with the present experimental results. Point LZ is the result of measurement 
by Schneider and Burton,“ made at 5 Mc/s by observing the penetrating 


13) Weber, Goeder: Phys. Rev., 61 (1942), 94. No description on frequencies. 
14) Schneider, Burton: J. Appl. Phys., 20 (1948), 48. There is no description on 
temperature, Presumably, the measurements were made at a room temperature. 
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power of supersonic wave through a plate at various incident angles. Dotted 
curves represent the recalculated values from data which have been presented 
by Alexandrov and Lazurkin ™ in terms of deformation amount of a test 
piece with definite dimension due to the oscillating external force with 
constant amplitude. Figures 1, 10, 100 and 1000 on the dotted curves show 
the frequencies per minute. The calculations were made by assuming that the 
amount of 4% of the strain corresponds to E=1.2x10’ dyne/cm2.46 However, 
whatever the assumption may be, there was an inevitable discrepancy between 
the results given by Alexandrov and Lazurkin and those given by Robinson, 
Ruggy and Slantz, the cause of which still remains unsolved. 

Though Fig. 7 seems to be a random collection of the fragmental results 
given in literatures, it may still serve to obtain some information on tempera- 
ture and frequency dependences of Young’s moduli. 

An inspection of Fig. 7 reveals the facts that: 

(1) There is a certain temperature range where a marked dependence 
of Young’s moduli on temperature as well as on frequency is observed. For 
example, at 100°C the dynamic Young’s modulus (at 100 cycles) shows the 
value of about one hundred times greater than that of the apparent instan- 
taneous Young’s modulus. These behaviours may be attributed to the elastic 
dispersion. 

(2) Outside of the said temperature range, temperature and frequency 
dependence of Young’s moduli seem to be rather slight. The steady value of 
Young’s modulus in the lower temperature region is about one thousand 
times greater than that in the higher. These mechanical behaviours of 
polymethyl methacrylate may be expected naturally by treating them as 
relaxational phenomena. 

(3) In the simplest model (mentioned in Part 1 of Section III) repre- 
senting the main features of the mechanical behaviour of high polymer, the 
total stress which is in equilibrium with the external force, is divided into 
two partial stress, one fi (purely elastic) being distributed on the primary 
bond and the other fs (relaxational) on the secondary bond. Let £\ and Ez 
be Young’s modulus of the primary and the secondary bond respectively.’ 
The total elongation under constant load is characterized by &, while the 
instantaneous elongation by E£,+:. The reasonable value of £; estimated 
from Fig. 7 is of the order of 10’ dyne/cm?.% At present, the knowledge 


15) Alexandrov, Lazurkin: Acta physicochimica, 12 (1940), 547. 

16) In order to bring the dotted curve 7 to the group of curves B and C as close as: 
possible, the said correspondence is made tentatively. 

17) It may be considered that the temperature dependence of those is rather slight. 

18) It may be remembered that the value is of the same order of that of raw rubber 


at a room temperature. 
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concerning the value of E(=E,+E,) is limited to the fact that it should be 
greater than the value of 6x1010 dyne/cm?, and its exact value is left 
undetermined on account of the lack of measurements at considerably high 
frequencies and at low temperatures. In the simplest model mentioned above, 
it is considered that the value of the partial stress attributed to the secondary 
bond would approach asymptotically to zero with time. The propriety of 
this speculation also still remains undetermined owing to the lack of facilities 
in making measurements under extremely long-time load application. To 
determine the propriety by measurements of smaller time scale, the measure- 
ments should be made at higher temperatures whereby distinct irreversible 
flow would be observed and difficulties of its separation from the total 
deformation would not be overcome.2 


V. Correlation between Transient and Steady-State Behaviours 


Part 1. 

In Part 1 the correlations between the transient and the steady-state 
behaviours are treated by a method which is suggested by that of Gross. 
Mathematically, the latter problem is completely analogous to the steady- 
state response of an electrical network to an alternating voltage source. 

Assume that the total deformation y(t under constant load Py is expressed 
by the following equation : 


> の = 1+ (4) Po (27) 


where 
£;: instantaneous Young’s modulus, 
F: form factor, 
W(t): creep function. 


The deformation y(t) due to the application of variable load P(é) is given 
by Gross with the aid of the principle of superposition as follows: 


tise lei +{ P@ he Cha dr| (28) 


where p(t) is the rate of creep. 
When the material is subjected to the sinusoidal stress, i.e., 


19) In practice, it is difficult to distinguish flow from delayed elasticity. It is not 


reasonable to regard the deformation as flow because of its apparent linear behaviour 
with time. 


20) loc. cit. 
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the steady-state admittance G(iw) is expressed by 
: F ; 
Gi =, | 1+AGw)—iB ©) | (29) 
t 


where 
y 4)=G Go)-P(£) 


and 


A (w)=|"cos (wT)- P(r) dt, 

0 

¥ (30) 
B =| a 

0 


Similarly, the steady-state solution of stress for alternating deformation 
is given as follows, i.e., assuming that the stress P(t) at fixed deformation 
yo is expressed by 


E, 


Pas 1-¥ © bo (27') 


where 
W(t): relaxation function, 


one obtains with the aid of the principle of superposition 


P= [2w-| yr) Padr | (28) 
where 
P (N= a : rate of relaxation. 


Putting y()=y-e'* and PW の =Z(iw) yD, the steady-state impedance 
Z(iw) is given by 


Zoe ze ses (oy Go@) | (29') 


where 


Aw)=| cos DSI 

(30/) 
B(o)={ sin (wT)-p (tr) dr. 

0 
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Expressing the complex elastic modulus by 
E* (—io)=E!+iE" (31) 
where 


E': dynamic modulus, 
E": elastic loss, 
one gets the following relations between modulus and impedance: 


F'=E,(1-A), EV=E)B. (32) 


The relation between impedance and admittance functions is given by 


G(iw)-Ziw)=1 (33) 
or by algebraic form: 
~ AC+A)+B2 5 B 
= | erent 34 
A (1+A)?+B? ’ (1+ A)? + B? (34) 


These correlations between the transient and the steady-state behaviours 
obtained by Gross lead to the deduction of the latter from the measurements 
on the former, i.e., when g(?t) is obtained from the experimental results on 
creep, A and B are calculated from Eq. (30) whereby A and B are computed 
from Eq. (34). The discussion below refers to the frequency dependence of 
Young’s modulus of ‘a substance, the creep of which is described by the 
following equation, a modified form of Eq. (1’), 


y(f)=at+b-InW-.). (35) 
Eq. (35) can be rewritten in the form: 


9 D=y (O)] 14+ on | (36) 


where 
y(0)=a+b- Inv. 
By comparing Eq. (36) with Eq. (27), one obtains 
(の = 0) 2 ge (37) 


Therefore, 


b 1 
t = = e => 
ptt) yO) (102%, t+o0 = t4+-0 (38) 
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The value of 8 is obtained as the slope of the straight line showing the 
relation between y and Int, while the value of y(0) is left unknown so long 
as the measurements of creep is concerned. However, suppose that the value 
of the instantaneous Young’s modulus E;(++E,+E2) assumes the smallest 
value (61019 dyne/cm?),2) the estimated value of y(0) is calculated from 
the following formulas: 

B 


ope: vs 
(0) = E; Pop CPS ABT (39) 


where 


I: moment of cross section of the beam. 
Zl: effective length of the beam. 
For example, from the experimental data that the creep curve of a beam 
with 7=25 cm and 1=(0.477)3x1.50/12 exhibits a slope of b=2.6x10-2 under 
constant load P=500x980 dyne, one obtains 


F 2.4104 
= — =e =0.20 
yO) E Py &x 1010 x 500 x 980 
Hence, 
Ge) 4345 = 00050 
y(0) 


Putting y(t)=C/(t+2), Aw) and BCw) are calculated from Eq. (30) as 
follows: 


Acay=c {Sea di=C\—cos wi: Cod) —sinad si Coa} 
tenis (40) 
Bo)=C | SE dt=C{—cos od «si (od) +sinod - Ci Coa) 
o Ott 

The results of numerical calculations of A/C and B/C are plotted in 
Fig. 8. Combining Eq. (40) and Eq. (34) A is obtained whereby the value 
of E'/E; (ratio of the dynamic Young’s modulus to the instantaneous one) is 
calculated from Eq. (32). The calculated values are plotted in Fig. 8, putting 
O=Land C=0.1. 

Fig. 8 gives a suggestion to determine, (1) the propriety of the specu- 
lation that, in order to represent the creep behaviour in.shorter time region, 
the factor Inf can be replaced by ln(@+ め , (2) the reasonable value of #, by 
comparing the experimental results obtained from the measurements of 


21) cf. p. 34. 
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Young’s moduli over a wide frequency range with the curve of 1— A4A(=E'/E;). 
However, there is so far no experimental result to refer except those given 
in the present paper (4.631010 dyne/cm? at 500 cycles, 4.4810! dyne/cm?2 
at 200 cycles) and the value of Schneider and Burton (6.6410! dyne/cm? 
at 5 Mc/s). Putting F;=6~7x10" dyne/cm2, a tentative attempt on the 
estimation of # from the above experimental data gave the values of 10-7~ 
10-8 sec, but in those cases the deviation from the curve was never less 
than 10%. 


Part 2. ; 

According to Dunell and Tobolsky,@ Kuhn had shown how, in case of 
a material for which the creep curve is linear with logarithmic time, one 
might predict the dynamic losses in free vibrations. They have pointed out 
22) Dunell, Tobolsky: J. Chem. Phys., 17 (1949), 1001, letter. No description on the 
deduction in their letter. : 

23) Kuhn, Kiinzle, Preissmann: Helv, Chim. Acta, 30 (1947), 307; 464; 839. Up to 


the present (1950) no further information is obtained except brief abstract on Chem. Abst., 
41 (1947), 3652; 4945; 5743. 
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that the rectangular distribution function of relaxation times, in terms of 
which Kuhn described an approximate mechanical system, had previously 
been extensively investigated by Becker@) for ferromagnetic hysteric be- 
haviour, and also that if stress relaxation is linear when plotted against 
logarithmic time, by using the same function, the following simpler corre- 
lation between stress relaxation and dynamic properties is obtained: 


TE 


Hf ieee Bo» 
ee 2 eo log 0 


x (slope of relaxation curve plotted vs. logy) time) (41) 


where 


1/bXwXVa, 


ayn: dynamic viscosity coefficient, 
€9: fixed strain at which the stress relaxation measurement is 
carried out, 
@: angular frequency, 
a, b: constants which are, respectively, smaller than and greater 
than the smallest and the largest recorded times. 

According to Weiss, the ferromagnetic substance contains a large number 
of small domains having intrincic value of magnetic intensity. In the absence 
of an external field, the resultant magnetization of the entire crystal is zero 
owing to the random orientation of domains. The magnetization of the 
crystal due to the external magnetic field is attributed to the instantaneous 
and retarded orientations of these domains, and it is considered that the 
principle of superposition is valid for retarded behaviour characterized by the 
distribution function of retardation times. Moreover, experiments on magnetic 
behaviours have revealed that the magnetic induction curve Bt) is linear 
with logarithmic time over a broad time range and the loss factor is inde- 
pendent of frequency over a broad frequency range. Studies on the mechanical 
properties of polymeric substances have also revealed that the transient 
behaviours are linear with logarithmic time and that the elastic loss is rela- 
tively independent of frequency over a considerable frequency range.‘ 

Thus, the magnetic behaviours of ferromagnetic substances seem to be 
completely analogous to the mechanical behaviours of high polymers, and yet 
little attention has been paid to the analogy. 

Owing to the lack of facilities in obtaining information concerning details 
of the theories of Kuhn as well as Dunell and Tobolsky, an attempt on the 
derivation of the said correlation is made on the basis of Gross’ theory, by 


24) Becker, Doring: Ferromagnetismus (J. Springer, Berlin), (1938), 103. 
25) A. Gemant, W. Jackson: Phil. Mag., 23 (1947), 950: Dillon, Prettyman, Hall: 
J. Appl. Phys., 15 (1944), 309: cf. letter of Dunell and Tobolsky. (loc. cit.) 
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using a frequency function of the rectangular form and also the function (24). 


Case 1. Stress relaxation at fixed strain. 


In Gross’ theory in which the creep is treated by a continuous system of 
Voigt’s model while the relaxation by that of Maxwell’s model, the mathe- 
matical expressions describing the relations between the creep and the 
dynamical properties are quite similar to those between the stress relaxation 
and the dynamic loss. Hence, assuming that the principle of superposition is 
valid, the relaxation function Wt) is expressed in terms of the distribution 
function F(7r) of relaxation times 7 as follows: 


y= 1— Fe ar | as 


where 


B: ratio of partial elastic stress to the final value of stress 


relaxation. : 
Therefore, the rate of relaxation is expressed in terms of relaxation 


frequency function p(s) as follows: 


dy 


= = 


|» (s)e-*s ds (20/) 
0 


where 


s=1/t: relaxation frequency, 
B For) dr=—p(s) &. 
From Eq. (20’) and Eq. (30’), B(w) is expressed as follows: 
ee eg all ods 
Bow)=|" pts) SS. (42) 


Suppose that the distribution function of relaxation times is of the form: 


1 i 
r=} A RT TD 


0 » (1 > 7 TT) 


(43) 


where 1/In(72/7;) is a nomalization factor. 
The corresponding relaxation frequency function p(s) is given by a 
rectangular form: 
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ee (Sh) 
as | In (42/71) T2 Ty 
B(s)= (44) 
labo 
T2 T1 
Combining Eq. (44) with Eq. (42), within the range of 
Th << i ぐ T2 
@ 
the function B(w) is expressed by 


Blo)= aa af f 2 = ods XO 8 as (45) 


In ( T2/T1) 1/72 @2+s2 In (T2/T 1) 2 


From Eqs. (45) and (32), the elastic loss E” is obtained as follows: 
B= Bos RB 0 ( 
SE iS 


Eq. (46) shows that when the relaxation frequenecy function is of a 
rectangular form, the elastic loss is independent of frequency within the range 
I 
of 1 oe. ぐ To. 
Next, combining Eq. (19’) with Eq. (43), within the range of 
TE Kt X€ T2 


one gets the following relaxation function for the rectangular relaxation 
frequency function: 


: 1 a 
Fi) Al 1+ Wer ee ae \ an 


where 0.577...... is Euler’s epee 
Substituting Eq. (47) into Eq. (27’), the total stress P(t) under constant 
deformation is expressed by 


he 2a 0.577--: In T2 
roastsRt I+ In (2/71) ‘es (72/71) H 


EO NR a eye (48) 


If the stress ere is found to be linear with logarithmic time, i. e., 
“P= =a’ —! Inz (49) 


one may make the correspondence between Eq. (48) and Eq. (49) in which 
the constants a’ and 8’ in Eq. (49) are expressed as follows: 
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hie ¢ 0577) « la (50) 
a’ =o gh tn GD GD | 


(51) 
las. Fr a 
8 =Yo i8 In (72/71) 
Eq. (51) can be rewritten in the form of 
B gi 
E; = —= (52) 
"In (72/71) ~ Yo 


From Eq. (52) and Eq. (46), the relation between the elastic loss E” and 
the slope of relaxation curve 8’ plotted against Int is obtained as follows: 


E" RS ees. (53) 
yo 2 


When the relaxation curve is plotted against logt, the following equation is 
obtained. 


jf Er sea (53!) 
2 yo logio 


Case 2. Creep under constant stress. 


If the measurements on creep behaviour are shown by Eq. (1’), similar 
mathematical treatment gives the following relation between the imaginary 


component of compliance B/E; and the slope 6 of the creep curve plotted 
against logt: 


al T 
7 B= 3P og (54)26’ 
Thus, when the transient behaviours are linear with logarithmic time 
over a considerable time range and the principle of superposition is valid, 
the simple correlations between the transient and the steady-state behaviours 
are derived so long as the frequency functions are of rectangular form. If 
the form of frequency function is flat over a considerable time range, it may 
be replaced by an appropriate rectangular function without significant error. 
It seems therefore most reasonable to conclude that the correlation mentioned 
above still remains to be approximately valid over the corresponding range. 
In part 2 of Section III, a distribution function of retardation times is 
given, in terms of which the linear relation between the creep under constant 
load and logarithmic time can be described within the ordinary time range. 
If the contribution of ¢ is neglected, one gets from Eq. (40) 


26) No literature concerning the frequency independence of imaginary component of 
compliance has been given. 
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lim. B(o)=C -—. 
o>) 2 


Therefore, within the range of ? « — «¢, the imaginary component of 


compliance is expressed by 


oa be ng 30s 
E; may ha ey ea Sie: y(O)logio att 
sts wb 
a 2P (0 )logi0 © Gy 


In view of the fact that Eq. (55) is of the same form as Eq. (54), it 
seems most reasonable that the two constants Y and ¢ which chracterize the 
function (24), correspond respectively to retardation times 7; and 72 in the 
rectangular function. When the contribution of ¢ is taken into account, the 
function B(w) is of the form: 


B(@)=C{—cos @ 3+ si(@v)+sinw? - Ci (oY) 


+cos wf - si (wf )—sin wf - Ci(wf)} 
where 
1T 
li i Ses, 4 
ut s1(x ) 2 


J limsi(z)=0, 


Xo 


lim sinx- Ci(x)=0, 
x っ 0 


him Ci(4)=0. 


x00 


Assuming ? < = « €, an approximate form of B(w)XC. i is obtained, 


and Eq. (55) still remains to be approximately valid. 
By a quite similar treatment on the'case of linear stress relaxation (Eq. 


he eS iy i 5 
(49)), the following results are obtained within the range of Y <ー—« és 


Y and ¢ being replaced by the corrsponding constants # and ¢, respectively. 


(56) 


Rk: 
E"=E; Bio) =E; Cie aes! (loge). Ps 


dae | tet 
~ 2y(0) logio 


where 


an 


1 
P(0) O36 


CR 
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Eq. (56) is identical with Eq. (53’). 

Comparison between the loss factors, predicted from Eq. (53’) (stress 
relaxation) and from observed values, was made on rubber-like substances 
and textile fibers by Dunell and Tobolsky. In most of the cases, the losses 
predicted from stress relaxation data were less than the observed dynamic 
losses. 

In order to make the similar comparison between the imaginary com- 
ponents of compliance predicted from Eq. (54) (creep) and from the observed 
value, it is required that the complex elastic modulus £;(1—A)+zE;B obtained 
by dynamic experiment is tranformed into B/E; by the following transfor- 
mation formula which is valid for any type of distribution function, i.e., 


Ria. CE; Bl 
E; fy CE; (1—A)? a LCE;B J 1 


In ordinary case, the second term in denominator can by neglected with- 
out significant error. 

Attention should be paid to the fact, for example, that when the trans- 
formation is made on the rectangular function p(s), the obtained function 
p(s) is generally of less flat form. Hence, it may by possible that Blo) is 
independent of frequency while B(w) is less independent of frequency, and 
that the time range, where the relaxation is linear with logarithmic time, 
is fairly broader than that where the creep behaves linearly. (Transformation 
problem is not treated here.) 

The verification of reliability of Eq. (54) by the use of data in literatures 
on polystyrene, polyvinylchloride, polymethyl methacrylate and ebonite has 
not proved successful. Since the data in literatures are those obtained for 
different purposes, it is thought advisable to carry out fitting measurements 
to answer the present purpose. By using the experimental data obtained on 
one and the same specimen, the comparison is made as shown in Table 1. 

From Table 1, it may be said that, except in case of polystyrene, there 
are still some discrepancies between the observed and predicted values. 
These discrepancies would arise from the difference of time scales in the 
measurements. i.e., in the case of creep, it took about one minute at least 
between loading and the first reading, while the measurements of the steady- 
state behaviour were made at frequencies of one half cycle or higher. 

In general, particularly in case of rubber, the predicted values are a little 
smaller than the observed, the results being similar to those of Dunell and 
Tobolsky mentioned above. The’ elastic loss G” of rubber increases with 
frequency even in lower frequency range,” while the rigidity G’ maintains 
its value whereby the value of G”/(G”"+G"’2) increases with frequency (G’ » 


27) Dillon, Prettyman, Hall: /. Appl. Phys., 15 (1944), 309. 


Viscoelastic Properties of Polymeric Substances 45 


Table 1. Comparison between the imaginary components of compliance 
predicted from Eq. (54) and from observed values. 


| Imaginary part of 
Material | Deformation | Temper- Frequency | compliance 
| | ature °C | —— 3 
| | Observed Predicted 
sees aa Se eens | ee Pa : ee 
| | 430 1.1 x 10-12 
Polymethy! | Bending = Trae Pike rages Pe A Cal 
methacrylate © (== ; | ーー 
| Torsion 20 iW 8.9x10-12 ¢ | 4.7x 10-12 a 
Polystyrene | Torsion 25 1/2 Te 410512) 6 | 1.5x 10-12 a 
Rubber (filer) | Torsion 30 1/2 7.9x10-9 ¢ | 2.4x10-9 a 
Rubber string | | iy. 
(Commercial | Elongation 20 1.96 | 8.3x10-9 d|2.7x10-9 e 
product) | | 
Phenolite | Torsion 27 1/2 | 1.0x10-M ¢ | 4.4x10-1 a 


a: cf. Section II. 
b= cf Section IV: 

: Measurement was made by means of multiple pendulum oscillator; S. Iwayanagi, 
T. Hideshima (cooperators in the present investigation): Presented before the 
47-th Scientific Lecture Meeting of Scientific Research Insititute, Tokyo. (1951). 

: Measurement was made by means of inverted pendulum; unpublished. (1952) 

é: unpublished. (1952) 


Attention should be paid to the fact that in general, particularly in case of rubber, 
the deformation in the case of creep is apt to be greater than that of vibration. The 
measurements on rubber were made under dead deformation. 


Notes on materials: 

Rubber: Manufactured by the Research Section of Sumitomo Denko & Co., Ltd.; filler 
15%; hardness 44; test piece was made ina cylindrical rod, both ends being cemented 
on brass disks. 

Phenolite: Supplied by Kawada, Tsuji Laboratory, Scientific Research Institute, Tokyo; 
condensed by ammonia catalyzer and sufficiently hardened, 


G"). Hence, a better agreement between the predicted and the observed 
values may be expected when the calculation is made by using the value of 
G"’ corresponding to lower frequencies. 

In case of polymethylate, the predicted value is greater than the observed 
in bending, while the former is smaller than the latter in torsion. In view 
of the fact that polymethyl methacrylate exhibits the secondary anomalous 
dispersion at a room temperature,™ it seems reasonable to attribute this 
inverse behaviour to the making of steady-state measurements at different 
frequencies. ‘ 


28) S. Iwayanagi, T. Hideshima: Presented before the Division of Polymer Physics 
at the Meeting of the Physical Society of Japan, Tokyo. April 1951. 
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In case of phenolite, the linear behaviour of creep with logarithmic time 
was somewhat irregular. 
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Classification of Snow Flakes and their Structures 


By 


Choji MAGONO and Hachiro OGUCHI* 


Abstruct 


Many snow flakes were photographed between 1951 and 1954 in Japan, 
and they were classified, according to the purposes by three methods, their 
metamorphic stages, shapes in falling state, and crystal habits. 

It was observed that the snow flakes change their crystal states in the 
process similar to that of the fallen snow. Considering aerodynamically the 
shape and velocity of snow flakes in falling state, a clue for presuming the 
uniformity of crystal habit in snow flakes was acquired. In the photographs 
of snow flakes obtained hitherto, it was found that in all types of snow 
crystals is formed one snow flake with the same type only, while another 
with almost all combinations between the different types is formed. All the 
ways of contact geometrically possible, namely, point, line, intertwined, and 
irregular contact were found in the snow flakes. The top of branch of 
dendritic crystals plays a leading role in the adhesion mechanism between 
the crystals. 

From those results, it was presumed that the large snow flakes might be 
composed of snow crystals of dendritic type. 


§1. Introduction 


Recently, some attentions have been payed to the snow flakes from the 
points of view that the snow flakes are considered to be the origin of large 
raindrops, and to correspond to the hydrometeors in the bright band of 
radar echo. 

Seligman? refered briefly to the manner of clinging between snow crystals 
in his “Snow Structure and Ski Field”. Many photographs of snow flakes 
were lately recorded in Nakaya’s “Snow Crystals”2) which was very useful 
to our work. Mason» described that the snow flakes originate from horizontal 
collisions between snow crystals during fall. One of the authors’? showed by 
calculations that a large snow flake could be formed by its vertical collisions 
with snow crystals due to the difference between their falling velocities. 

Although theoretical considerations have not been yet made, it appears 
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that the snow flakes are formed near the freezing level, and the snow crystals 
of dendritic type are easy to form the large snow flake. In order to investi- 
gate the snow flake systematically, it may be necessary to act on the results 
of numerous observations. With this aid, we have gathered as many photo- 
graphs as possible in recent four seasons in Northern Japan, and arranged 
them on the basis of their structures for the next process of investigations in 
this field. 


§2. Classifications of snow flakes 


According to the purposes, snow flakes were classified by three methods, 
that is, their metamorphic stages, shapes in falling states, and the crystal 
type of snow crystals which occupy the majority of the snow flakes. 

a) Classification by metamorphic stages. 

Fresh snow flakes suffer various changes in their falling states, to say, 
rimed by supercooled cloud particles, granulated by sublimation or refreezing, 
‘and wet by partially melting, since then melted into raindrops at a warm 
climate. The changes correspond to those of fallen snow which Saito® classi- 
fied by their metamorphic stages. Our classifications of snow flakes are, there- 
fore, more’or less similar to that of Saito, as stated later. The snow flakes 
at each stages are shown in PI. I. 

’ ‘Fresh snow flakes (Photos. 1~3, Pl. I). Inacold climate, one may observe 
snow flakes whose snow crystals do not suffer remarkable changes yet and 
hold fine structures in their branches or connecting parts. This is the initial 
stage of the snow flakes. We call them in this stage “Fresh snow flakes”, 
neglecting a few cloud particles attached on the snow crystals, because it is 
rare near the earth’s surface that snow crystals carry no cloud particles. 

Rimed snow flakes (Photos. 4~6, Pl. 1). Such snow flakes that are covered 
by numerous cloud particles, we call “Rimed snow flakes”, although it is 
unclear that after forming the snow flakes the particles have attached to or 
initially rimed crystals have clinged mutually. 

Granulated snow flakes (Photos. 7~9, Pl. 1). When a snow flake falls 
into an atmospheric layer of relatively low equilibrium vapor pressure, rapid 
sublimation takes place, and the snow flake loses the fine structures in its 
snow crystals. At the lower layer, individual crystals of the snow flake 
undergo also melting and refreezing partially. Owing to those effects, the 
crystals tend to be granulated. This is the second stage of the snow flake. 

Wet snow flakes (Photos. 10~12, Pl. I). When a snow flake is considerably 
wet, it shrinks undoubtedly but does not become always one block, In the 
case of Photo. 11, the snow flake has been separated to three parts. Consider- 
ing that the snow flakes are constantly subjected to the resistance from air in 
the vicinity during fall due to their large falling velocities and .complicated 
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shapes, it is supposed that snow flakes would be disrupted into a few raindrops 
before melting perfectly. 

Melted snow flakes (Photos. 13~15, Pl. I). This is the final stage of snow 
flakes. In sleety weather, a piece of ice unmelted yet is detected within a 
raindrop, although this ice piece is often overlooked. Photos. 13-15 show 
such pieces contained in raindrops by which some splashes were raised, as 
seen in those figures obtained by photographic paper method.® 

b) Classification by shapes of snow flakes during fall. 

It is our opinion that the snow flake in falling state should be more 
watched, because the name “snow flake” is given to the snow during fall. 
From this point of view, we classified the snow flakes by their shapes while 
falling, as follows. 

Horizontal type (Photos. 16 and 17, Pl. II). It is aerodynamically stable 
that a body of plate type or a pole type with a uniform density, falls keeping 
its plane or axis horizontally. Taking an analogy with this phenomena, we 
presume that such a snow flake that keeps its plane or axis horizontally 
during fall, may be composed of uniform snow crystals, in other words, have 
a uniform density. For examples, Photo. 16 shows a snow flake of horizontal 
type formed of uniform dendrites, and Photo. 17 represents that of two stellar 
crystals. It is considered that the snow flake of this type has a larger 
probability of touching mutually than those of other types, because the former 
occupies a wider cross-section than the later and falls slowly, that is, floats 


for a long time in the air. 
Vertical type (Photos. 20 and 21, Pl. II). If one end of a snow flake is 


extremely dense, the snow flake falls keeping its heavier part downward and 
fluttering its lighter one. Photos. 20 and 21 are the examples of this type. 
These snow flakes are composed of a graupel and a stellar crystal which are 
considered to have clinged mutually owing to the difference between their fall 
velocities. From observations,” it was found that the snow flake of vertical 
type has larger fall velocity than that of horizontal type. 

Inverted cone type (Photos. 18 and 19, Pl. Il). The snow flake of this 
type is situated between those of horizontal and vertical type. During fall, 
the relatively heavier part of the snow flake occupies the downward of the 
cone, and an end of the lighter one plays a rudder with a constant angle. 
At the result of the rudder, the snow flake of this type falls rotating spirally. 
One may see the rotating motion and the compact part at the downward of 
inverted cone in Photos, 18 and 19. 

c) Classification by crystal types. 

In our observations, it was found that all kinds aie snow crystals except 
sleets, form snow flakes so we classified also the snow flakes by the classi- 
fication of crystal type of component snow crystals for practical purposes 
by the International Commission of Snow and Ice. After this, we shall repre- 
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sent the type of a snow flake according to that of snow crystals which occupy 
the majority of the snow flake. Sleets (refrozen ice pellets) are scarecely ob- 
served in our country, and it is not likely that they form snow flakes. Observed 
snow flakes are tabulated in Pls. III and IV, and the explanations of each 
photographs are given in the opposite pages of the corresponding plates. 

The base of the snow crystal of columnar type shown in Photos. 22, 23, 
and 24 have grown somewhat horizontally, namely, have approached to capped 
columns. We think that pure columns do not form recklessly snow flakes, 
although crystals of bullet type make cling one another at their tops as 
descrived in “Snow Crystal,” p. 426. 

The classifications described above are tabulated in Table 1. 


Table 1. 
Classifications of snow flakes. 
eee Term Code Photo. no. | Pl. no. | 
fresh snow flake CK 1-3 | I | 
rimed ” ” CK) 4-6 | I | 
Metamorphic stage granulated » " (O) 7-9 . I | 
wet ” ” ™® | 10-12 | I 
melted 4 | 13-15 I 
; horizontal type | ° 16,17 | 
oles falling inverted cone # Vv 18,19 | II 
vertical ” 0 20, 21 II 
columns | 2208 | III | 
plates 25, 26 Ill | 
stellar type 28-30 III 
spatial dendrites 31,32 Ut 
Crystal type capped columns 34, 35 IV 
needles 37, 38 IV 
; graupel 40, 41 IV | 
| irregular crystals 43, 44 IV | 


§3. Combination between snow crystals of different types 


Snow crystals of different types are often observed at a same time, and 
those of each types must have various falling velocities respectively. If. the 
snow flake is formed by the collision with the snow crystals owing to the dif- 
ference between their fall velocities, it is probable that a snow flake collides 
more frequently with the crystals of other type than with those of same one. , 
And it is suggested that in a snow flake, a component snow crystal of a larger 
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fall velocity might grow at a higher atmospheric layer than the layer at 
which the other crystal of smaller velocity. In the fact, the snow flakes com- 
posed of various combinations between different types were observed as shown 
in the last columns of Pls. II and IV. For the simplicity, we selected only 
the combinations between two different types which are graphically shown in 
Fig. 1. The thickness of the lines which combine two types, represents the 
frequency of occurance determined by our photographic observations. There- 
fore, the frequency does not necessary to show the grade of the inclination of 
clinging between the types, but perhaps mean the frequency of occurrance of 
the snow crystals. This conception is supported by the fact that the frequency 
of snow flakes agrees with that of individual snow crystals. 


Types of snow crystals  ! 
composing snow flakes Combinations 


Columns 


Plates 


Stellar type (sector form) 


Spatial dentrites (radiating 
type) 


Capped columns 


Needles 


Rimed crystals (thick plates, 
graupel-like crystals) 


Graupel 


eo IL OX © 


Irregular crystals 


Fig. 1 
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The meteorological conditions for the growth of various crystal types are 
shown in Table 2. It is natually expected that snow crystals which grew at 
similar air conditions may form a snow flake, but it is remarkable that as 
shown in Photos. 20 and 42, a stellar crystal clings to a graupel whose 
conditions for growth is considered to be far different from that of stellar 
type. Such a snow flake might originate from the difference of fall velocity 
between the two crystals. One sees in Fig. 1 that the snow crystals of 
irregular type combine with almost all types of crystals. This may be caused 
by the following facts; it is often difficult to discriminate the crystal type, 
therefore, the unclear crystals tend to be regarded as of irregular type, and 
the crystals of this type fall more frequently than is expected generally. 


Table 2. 
Meteorological conditions necessary for growth and falling velocities 
of snow crystals of various types. 


Crystal type Air tompsoture Vapor supply ~ | ey | 
| Needles 1-4 : 0. 3-0. 6 
Irregular needles 4-7 . 
Columns 9-12, 17-21 little | 1 | 
Plates, sectors 11-14, 17-21 little 
Dendrites (stellar, - . 
spatial) 14-17 much 0.3-0.6 | 
| | 
| Rimed crystals : : | 
(thick plate) 9-14, 17-19 with cloud particles 0.6-1.8 
Graupellike crystals with cloud particles | 
Graupel relatively warm With cloud particles | 1-2.5 
Irregular crystals | 0.6 


§4 Manners of contact between snow crystals 


We classified geometrically the manners of contact or clinging between 
two snow crystals as follows, point contact, line contact, plane contact, 
intertwined contact, and irregular contact. These manners of contact were 
more finely classified as shown in Fig. 2 in which typical examples are 
sketched from corresponding photographs. Such a assemblage of snow crystals 
that are considered to have originated from the attachment of the nuclei i.e. 
twelve sided crystal or combination of bullet type gathering at heads each, 
were excluded from the snow flake. 

a) Point contact. 

This manner of contact is more classified to point-point, point-line, and 
point-plane contact. 

Point-point (Photos. 46~51, Pl. V). Almost all dendritic crystals cling 
to one another at the ends of thier branches. The examples are shown in 
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Photos. 46~51, and their explanations are given in the opposite page of the 
Pl. V. Frequently, crystals of stellar type cling spatially one another, thus, 
a very light snow flake can be formed from the stellar crystals, although 
the individual crystals are not so light. Some times. a large dendritic crystal 
of plane type carries many small crystals at the ends of its branches, and 
becomes often the center of a snow flake. See Photo. 50. Rarely two snow 
crystals cling mutually with such a manner that the direction of the branches 
of both crystals are parallel as seen in Photos. 48 and 51. 

Point-line (Photos. 39 and 35). Some rimed irregular crystals cling to 
fine needle crystals with point contacts. 

b) Line contact. 

Line-line (Photo. 25, Pl. III). Such crystals as plates or columns which 
grow at the condition of low super-saturation vapor pressure cling rarely at 
their edges. 

Line-plane (Photos. 22 and 23, Pl. III, Photo. 34, Pl. IV). Deformed columns 
are inclined to cling to one another with this manners of contact. 

c) Plane contact. 

Snow flakes composed of stellar crystals which were piled up each other, 
were not so frequently observed as is expected. The snow crystals shown in 
Photos. 28~30, Pl. III are apparently piled up one another, that is, touched 
with the plane contact, but they might contact mutually at the fine points 
on their surfaces. 

d) Intertwined contact. 

Snow crystals intertwined as sketched in Fig. 2 are descrived in “Snow 
Crystals,” p. 277. It is likely that these intertwined crystals gathered each 
other at a half way of their crystal growth, because it is not expected that 
these two stellar crystals intertwined by accident, leaving a especially strong 


adhesion power at the ends of the branches. This adhesion power will be " 


stated later. 

e) Irregular contact. 

Snow crystals included in large snow flakes contact irregularly by mutual 
mechanical pressure. 


§5. The adhesion power between snow crystals, 
forming a snow flake 


As descrived above, it is almost all at the ends of branches or tops of 
cones that snow crystals cling to one another. Because the ends of tops are 
situated at the most outer parts of snow crystals, it would be natural that 
snow crystals cling often at those parts if they contact mutually. 

But on the strong adhesion power at the fine parts, we think, some pro- 
blems are still left. It is surprising that a graupel of large fall velocity holds 
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still the touch with a fine dendrite during fall. 

When two fine points of branches contact with each other, the point 
would melt temporarily owing to the mechanical pressure or frictional heat. 
Then a strong bridge may be formed in the place by sublimation phenomena, 
as shown in Photo. 52, Pl. V, because the contacting part have very large 
curvature. 

In the case, snow crystals lose their fine structures and become round with 
the progress of sublimation, as Saito» and Kojima’ pointed out in fallen 
snow. Cloud particles which froze to the surface of a crystal would play 
also this fine connecting role. On the other hand, some snow flakes which 
contain rimed snow crystals are fragile when received on the ground. Photos. 
53 and 54, Pl. V show the fragments from such a fragile snow flake respec- 
tively. It is also considered that Coulomb’s force between snow crystals 
charged at opposite signs, has some effect on the adhesion power. As the 
electrostatic charge is localized to the corners of snow crystals, Coulomb’s 
force would be particularly strong at these places. The graupel which is 
formed of numerous cloud particles are considered to be charged negatively 
by Workman-Reynolds’ effect. While, it was found by many observers” 
that the charge on snow is more often positive than negative. Therefore, 
Coulomb’s force will probably play a main role in connecting the graupel 
and stellar crystal. 


§6. Large snow flakes 


Such a large snow flake that is composed of thousands crystals, we 
think, would have grown under following circumstances. 

i) The larger the concerned space density of snow crystals is, the more 
often the crystals will collide mutually. From the calculations by the one of 
the authors,” the volume of the snow flakes is proportional to the third 
power of the space density of the crystals in the air. The condition on which 
the rate of the crystal growth is largest, ranges from ~14 to -17°C, because 
the difference of saturation vapor pressure between super-cooled water and 
ice is most in the range. This temperature range agrees with that of the 
condition on which snow crystals grow to the dendritic form. 

ii) As the snow crystal or flake of dendritic type is very light, in other 
words, has a large volume per one crystal or flake, its chance for the mutual 
collion of this type may be more than of other types. 

iii) The falling velocity of dendritic snow flake is small, that is, the 
snow flake of this type floats in the air for a longer time than of other 
types. Therefore it has the largest frequency of collision with snow crystals 
in all types. The snow flakes of high fall speed can not grow largely, as 
descrived in § 2. 
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iv) It is considered that the snow flake begin to be form at about -10°C 
level. From this point of view, the snow crystals of dendritic type which 
grow at -15°C level, are favorable to a wide mutual collision layer, compared 
with those of other types which grow at the warmer i.e. lower level. 

v) As descrived previously, the dendritic snow crystals have many 
branches which have a strong adhesion power. 

vi) A calm weather and relatively warm climate near 0°C are necessary 
for the large snow flake growth, of course. 

Considering these reasons for snow flake’s growth, the large snow flakes 
must be of dendritic type, and of horizontal one, because the snow flakes of 
vertical or inverted cone-type can not grow largely due to their high fall 
speed. The large snow flakes shown in Photos. 55~57, Pl. V support our 
opinion. 
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Explanations of Plate II 


Side views of snow flakes while falling, photographed by successive 
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electric sparks. 


A snow flake of horizontal type, being fluttered, 

A snow flake of horizontal type, composed of two stellar crystals. 

A snow flake of iaverted cone type. 

A snow flake of inverted cone type. The lower part is more compact 
than the upper part. 

A snow flake of vertical type, composed of a Steller crystal and a 
graupel. 

A snow flake of vertical type. 
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Explanations of Plate III 


Types of snow crystals, composing snow flakes. 


Assemblage of columns; irregular contact. 

Three columns-plate ; line contact. 

Column-plate; point contact. 

Plate-plate ; line contact. 

Plate-plate ; point contact. 

Plate-column; point contact. 

Two dendritic crystals of plane type-stellar crystal; piled up contact. 
Sector-sector; point contact, (Photographed by Mr. K. Higuchi). 
Three rimed crystals of plane type; piled up contact. 

Spatial dendritites-stellar type; point contacts.” 

A part of a snow flake, composed of spatial (radiating) dendritites ; 
point contacts. 

Spatial dendrite-graupel; plane contact. 
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Photo. 28 


Spatial dendrites 
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Explanations of Plate IV 
Type of snow crystals composing snow flake. 


Capped columns; line contacts. 


x5.7 Capped columns-spatial dendrites; irregular contacts. 


x21 


Capped columns-plate; line contact. 


x 10.5 Needles; point contacts. 


x 22 
x12 


Needles; point contacts. 
Needles-rimed irregular crystals; point contacts. 


x5.2 . Graupel-graupel; point contact. 

x5.5 Graupel-graupel; point contacts. 

x9.5 Graupel-stellar crystal; point contacts. 
x10.5 Irregular crystals; irregular contacts. 
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Irregular crystals-stellar crystal irregular contacts. 
Irregular crystal-needle point contact. 
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Photo. 34 
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Explanations of Plate V 


Manners of contacts between snow crystals, 


Contact at an end of a branch, (photographed by Mr. K. Higuchi). 

A small crystal, clinging to a top of a branch, 

Parallel contact of branches between two crystal, clinging by one branch, 
(photographed by Mr. K. Higuchi). 

Three crystals, clinging spatially by two braches each. 

Dendritic crystal, carrying rimed crystals with it at ends of its branches. 
Parallel contacts of branches between two crystals, clinging with many 
ends of their braches, (photographed by Mr. K. Higuchi). 

A combining branch, thickened by sublimation, 

Fragments from a snow flake. 

Fragments from a snow flake. 

A large snow flake, 


Photo, 57, 58 1.3 A large snow flake in which many capped columns (H-shaped) are 


seen, 
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Study on Lichtenberg’s Figures by Means 
of Special Dry Plates. 


By 


Bunji ARAI 


Abstract 


The Lichtenberg’s figures were researched on various conditions of the: 
photographic dry plates, and these figures were obtained by applying a 
constant crest voltage of impulse. Special dry plates were adopted and 
researched as follows ;— 

It was observed that the change of the degree of emulsion ageing (or 
ripening) affected almost no variation on the positive and negative figures 
respectively. For the plates which added a desensitizer (including the nuclear 
track plates, the autoradiographic plates and the plates for electron beam), 
the existence of the desensitizer was dominant. There was observed great 
reductions of the size on both figures respectively. But the records of positive 
figures’ streamers were fairly observed on the plates for electron beam, and 
clear records of coronets were obtained. For the plates of various electro- 
conductivity (A), gained by adding KNO; to tke emulsion, there was observed 
small variation on the positive figures, and in the case of increased value of 
K, fine branched dendritic streamers appeared on the same figures. But, on 
the negative figures, remarkable enlargements of the figures were observed 
for the large value of K. At last, for the variations of silver halide (AgBr) 
content, there were observed enlargements of the size of the figures for the 
reduction of AgBr concentration on both figures, and vice versa. 


§ 1. Introduction 


Lichtenberg’s figures are varied by various factors. Their variations are 
caused by the electrical,” and atmospherical conditions?’ and the conditions 
of photographical emulsion etc. In upper conditions, the writer has not yet 
heard of the experimental studies made on the Lichtenberg’s figures by 
changing characteristics of the emulsion, except only C.E. Magnusson’s 
report?’ on the Lichtenberg’s figures examined on the various kinds of plates: 
available in the market in about 1930. This paper was described with the 
object of the preliminary experiments, aiming to raise the accuracy and 
extend the range of the measurements, when the Lichtenberg’s figures were 
- used as the methods of the measurements of an impulsive voltage. 
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§ 2. Special Dry Plates, Experimental Apparatus and Method 


Special dry plates were prepared as follows. 
(1) The plates which had the degree of emulsion ageing changed, they were 
a series of AP-155A (non ageing), AP-155B (medium ageing), AP-155C 
(max. ageing). 
(2) The max. ageing plates which added a desensitizer (AP-155D). 
(3) The nuclear track plates (ET-2E, their thickness of the emulsion layers 
were 50 /). 
(4) The autoradiographic plates (ET-2E, 15 gp thick). 
(5) The plates for electron beam (EMG 1). 
(6) The plates of various electro-conductivity (4), gained by adding KNO; 
to the emulsion. They were a series of AP-157A (kK: 1000x10-§2-cm~-}), 
AP-157B (KX: 3000 10-6.2-cm-1), AP-157C (AK: 9000x10-6.2-cm~-!), AP-157D 


ok RG high tension source 
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Fig. 1. Wiring diagram of impulse circuit. 
Cs: Large Leyden jar (0.0038 pF x2). 
Gi: Sphere spark gap (2.5 cm ¢ brass ball, fixed at 0.95 cm, 25°C, 760 mm Hg., 
30 kV—peak value). 

G2: Sphere spark gap for chopper of impulsive wave (2.0 cm ¢ brass ball, 

used at 0.35cm, 20°C, 760mm Hg., 13.8kV—crest voltage). 

R: Fixed high resistance (HOKUOHM—resistance ready-made as a radio 

parts, 1 MQ x10, 1/2 watt type). 
r: Liquid resistance for damper of impulse (saturated solution of CuSO: 
0.26 kQ, 20°C). 

Ly, Ly,......: Hand made Klydonograph camaras2) by the writer. 

Gy, G2 were illuminated with the ultra violet rays from a mercury vapour lamp, 
aiming to give the min. fluctuations of the spark discharges. The lead wires which 
from the terminals of G2 to the each cameras weré used same length one, and con- 
nected in parallel on each other. 


“to D.C. high tension source” leads to Gaiffe-Gallot et Pillon’s apparatus). 
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(K: 18000 x10-69.cm-1l), these numerical values of K were measured at the 
emulsion in the liquid state (40°C). 

The emulsion of AP was different from the emulsion of A1, but funda- 
mental properties were the same as those of A 1. 

(7) The Fuji Al plates in the market worked up by the writer. They were 
dipped (5 minutes in moving state) in the various concentrations (1%, 2%, 
5%, 10%) of NaCl solutions to give various electro-conductivities and were 
dried in desiccators in a dark room (20 hours). 

(8) The plates which were varied of silver halide (AgBr) content (a series 
of ET-2E, AgBr content was 71.0%, 35.5%, 23.7%, 17.7%, 14.2% and the re- 
mainded was gelatine). 

These special dry plates had been made in research laboratory of the 
Fuji Photo-Film Co. Ltd. and were favourably presented to the writer. 

For the comparisions were adopted the Fuji A1 plates. 

The apparatus is the impulse circuit hitherto used, which was supplied 
from the Gazffe-Gallot et Pillon’s D.C. high tension source”. See Fig. 1. 
The four or five kinds of special plates set up in hand made cameras» were 
connected parallel so that the same shape and the same crest voltage of 
impulse might be supplied at the same time in the positive or negative 
polarity. 

These plates were developed by the same developer (DK;) at the same 
time. 


(oan high tension source 


Fig. 2. The other connection diagram of impulse circuit. 

Connection to these cameras, to supply the same shape and the same crest 
voltage of impulse, to obtain the positive and negative figures on the same emul- 
sion surface at the same time. The notations are the same at that for Fig. 1. 
But, G2=0.78 cm, Rr=Rz2=1 MQ x5. 
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Also, the circuit shown in Fig. 2 was used owing to insufficient samples. 
With this wiring the positive and negative figures can be obtained on the 
same emulsion surface at the same time, the two cameras being in parallel 
connections. 

Through the experiments, the sphere spark-gap Gi(2.5 cm ¢ brass ball) 
was kept at 0.95cm (30 kV peak value, 25°C, 760mm Hg.), the impulse sup- 
plied by G, was cut off by the sphere spark-gap G2. G2(2.0cm $& brass ball) 
was used at 0.35 cm in the connection of Fig. 1 and 0.78 cm adopted in the con- 
nection of Fig. 2 (both case supplied 13.3 kV peak value, 20°C, 760mm Hg.). 
The reason why the peak value 13kV is supplied, is that this crest voltage 
gives a minimum mean error in the radius of the positive figures (45% )”. 


§3. Experimental Results and Discussion 


Photo. 1 (Pl. 1) shows the series of AP-155, which were obtained by the 
wiring of Fig. 1. It was observed that the few variations of the figures by 
the change of the degree of emulsion ageing were affected on the positive and 


negative figures respectively. But the effect of the desensitizer (Sdaure- . 


Rhodemaire) was dominant. There were observed great reductions of the size 
on both figures respectively. Especially, the AP-155D could not record the 
streamers of the positive figures, and obtained skeleton figures of the coronet 
only. Photo. 2 (Pl. 2) shows that there appeared the same results by the 
wiring of Fig. 2. For the nuclear track plates and the autoradiographic 
plates, there was observed the same tendency of the figures as shown in 
AP-155D in Photo. 1 and 2. But, in the nuclear track plates of the thick 
emulsion (50) obtained a smaller size, compared with the autoradiographic 
plates of thin emulsion (15). Photo. 3 (Pl. 3) shows the example of this 
result. In spite of the desensitizer existent, the plates for electron beam 
(EMG 1) recorded fairly positive figures streamers. Especially clear records of 
coronets were observed. Photo. 4 (Pl. 4) shows the results by the wiring of Fig. 
2. The left side shows the Fuji Al plates for comparisions at the same time. 
In these figures, the writer examined the inverse figures on the same plates. 

In the series of AP-157, the results of the positive figures are shown 
in Photo. 5 (Pl. 5). In the case of increased value of K, there were 
observed few variations. And fine branched dendritic streamers appeared 
on them. The results of the negative figures are shown in Photo. 6 (Pl. 6). 
In the case of increased K, the size of the figures remarkably increased and 
AP-157D flashed over. In this case, the writer thought that these facts 
moved easily the electrons of field emission from needle electrode for the 
existence of KNO; in the emulsion. Photo. 7 (Pl. 7) shows the other ex- 
amples of both figures. The right side column shows negative figures, aim- 
ing at comparission of Al plates which were cut out from one sheet. There 
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were observed no instrumental errors of each camera. In the case of AP-157C 
and AP-157D the crystals of KNO; were deposited in the back grounds and 
the black spots appeared in these parts when reprinted. 

Photo. 8 (Pl. 8) shows the examples of the Fuji Al plates treated with 
NaCl solutions by the wiring of Fig. 2. The left side figures shows the Fuji 
Al plates of no treatment for comparisions at the same time. In the case 
of increasing NaCl concentration, the streamers of the positive figures were 
branched in a fine dendritic (like from the examples of 5% NaCl concen- 
tration), and on the plates treated with 10% NaCl solutions were observed 
a few thin streamers. In this case small white stars appeared all over the 
surface when. reprinted. This fact would be a result from the deposited 
crystal of NaCl emitting fluorescence in a strong electrical field. In these 
plates the variations of the size of both figures was scarcely observed any 
variation. 

See Photo. 9 (PI. 9). This shows the result of the positive figures on the 
plates (ET-2E) for the variations of silver halide (AgBr) content. The 
numerical figures of % written in each figure show the amounts of silver 
halide content, and the remainder was gelatine. Fig. 10 (Pl. 10) shows the 
series of the negative figures. 

In the case of high concentration of AgBr was observed a decrease of the 
size of both figures, and vice versa. 

From the above observation, this fact would suggest that the extension of 
measurement range can be made by adjustment of silver halide content. 

This paper is a summary of the lecture made at the Special Meeting of 
the Applied Physics Department, held by the Physical Society of Japan, in 
Osaka, on Nov., 5. 1954. 
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Date and data for Photo. 1—10 


che Date Weather a seal Room Temp.(°C) RHA (%). Remarks 
1 | Oct. 24. 1954 | Cloudy 767.6 | 18.4 een eed ee 
2 | Oct. 21. 1954| Half fine \761. 1->760. 2 17.9 79 a 
3 | Apr.14.1950| Fine |765.6 >765.1 17.8 , 64-69, |ET-2E (50x) 
» | Sep. 1. 1951 | Cloudy 754.0 27.9 oF eS 
OR ER) ig 762.7 22.8 75. CS 
ah ae 科 A 762.5 26.3 -| - 65 - | -Al, EMGI 
5 | Oct, 25.1954] 769.5 18.0 75 | AP-157@ 
6 5 2 769.0 17.9 73 i 
FE : 767.2 18.3 2 AP AST ao 
Se get Cabem a 765.1 18.4 8 | AOS 
Bg | Sep. 23.1951; 4b. 767.0 23.5 | 88 . Re a x, 
on SORES sae 761.4 2559 iPS See A hae 
9 | Oct, 28. 1954| Half fine | 758.0 17.4 ey ET-2E ® 
10 st e 787.1 17, 6.doxghiorpis se cle aes 


* By wet & dry bulb hygrometer. 
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Photo. 1 


Postive fig. ~V=13.5kV, (Vy =13.3kKV)> Negative fig. 


AP-155 D 


| 
AP-155 C 
| 


AP-155 B 


AP-155 A 


Vy: Sparking peak voltage at 760 mm Hg, 20°C (Gp). 
V: Corrected Sparking peak voltage (G»). 
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Fuji ET-2E, 50u, V=4.2kV, (V,=4.1 kV). 
(Nuclear track: plate). 


Fuji Survey-plate. 


”» 


Fuji ET-2E.15u, V=4.0 kV, €V,=41kV). 
(Autoradiographic plate). 


Oriental 1200 plate, V=4.2kV, 
Vy, =4.1 kV). 


Vee WEKV)  IBOKV. 10.5 kV 9.0kV 7.4kV 


‘O) 


B. Arai: Photo. 4 PI. 4 


Fuji Al EMGI 


V=13.2 kV 
CV, =13.3%;, ) 


Fuji Al aren 


10.8kV 
(105°,, ) 


10.3,, 
105 


3 EMG1 


Pl. 5 


Ty tng 
u 0008T :C wv 
Ga 
‘s2nBiq ostysog (ANSEL = “A) ‘ANMNET =A : Ht ie : 
y-Wid-oyU 9-0OT X0001 + VLSI-dV 
WO Git FT ger Fee Wo ped a ae ; 
tr | | | ' i | 
fe) 
8 
a 
CLS Tak Vi 2 LST-dV q LS1-dV WPAN GE Se 


ie 


B. Ara 


© 
Ly Hay 


I 


a 
4 OCOsI -d / 
. » - Bie’ ‘ ° = 4 0006 0) 7] 
seinsy 9Al}2BoN CAN S'EL= “A) ‘ANMEL=A i b00e Seg 4 
t-W0-OYW 9-61 X000I :V LSI-dV 
a ZST-dv 
に =】 
8 9 ZSI-dV 
fe) 
Ay 


d LS1-dV 


B. Arai: 


V LSI-dV 


B. Arai: 


Photo. 7 


V=13.5 kV, CV, =13.3kV) 
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Pl. 10 
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Real Paths of Meteors Observed in August 1950-1954. 
by 
Sigeru KANDA 

Abstract 
The writer obtained real paths of 24 meteors from the observations made 
by several observers in August 1952, in the former paper. In this paper, 
from the visual observations of meteors made by more than 43 observers, 
including students of Yokohama National University and members of Japan 
Astronomical Study Association, in August 1950, 1952, 1953 and 1954, real 

paths of 79 meteors have been determined. 


先 に 筆者 は 「1952 年 8 AMMO WEO REM) LT SIP RR Lica, それ 
は 1952 年 8 AORMBMD 6 24 個 の 流星 の 実 経路 を 決定 し た も の で あつ た 。 本稿 
は その 研究 の 続き で あつ て , 1950, 1952, 1953 年 及び 1954 年 の 8 Atco Sei 
観測 か ら 決 定 し た 流星 の 実 経路 の 研究 で ある 。 

第 1 実は 実 経路 を 決定 し た 流星 の 観測 者 の 観測 地 並 に その 経緯 度 を 示し た も の で , 
EEO MUA SAU. 観測 地 は 宮城 , 新潟 , . 栃木 , 埼玉 , 東京 , 神奈 川 , 静岡 
の 都 県 に 芯 つ て いる 。 個人 観測 者 25 名 及び 8 団体 で ある が , 団体 観測 者 で も , 個人 
別に 観測 を 報告 され て いる も の も あり , それ 等 は 実 経路 の 決定 に 関係 の ある 氏名 の み 
WFAA UK, 観測 者 の 総数 は 43 名 以上 で ある 。 都合 に より 前 報 と 路 符 を 改め た も の 
が ある 。 

第 2 表 は 実 経路 を 決定 し た 流星 の 各 観 測 者 の 観測 で ある 。Wt. は 観測 の 重み で 1~ 
5 の 5 階級 こ による 。 Dur. kee, Mag. teocHf, Vel. は 速度 で 繧 か の 方 か ら 
vS, S, 18, M, rR, R, VRO 7 HMchS. Bt OK H YUH, WEA, RE 
ds, BCH. Biginning は 流星 の 出現 点 , Ending は 消 減 点 で , 赤 経 ¢ 赤 緯 8 
を 以 て 示し た 。 TRROSOKRLLRAL, 附 記 し た 数 字 は 療 の 残 つ た 秒 数 で ある 。 

実 経路 を 決定 し た 流星 の 総数 は 1950 年 32 個 , 1952 年 13 個 , 1953 年 27 個 , 
1954 年 7 個 , 計 79 個 で あり , 前 報 の 24 個 を 加え れ ば 総数 103 個 と な る 。 流量 の 
番号 は 1950 年 の も の は 5001~5032, 1952 年 の も の は 前 報 の 24 個 に 続い て 5225~ 
5237, 1953 年 の も の は 5301~5327, 1954 年 の も の は 5401~5407 と いう よう に , 一 
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見 し て 観測 の 年 が わか る よう に つけ た 。 

第 3 RIL 2 BOMMD ら 決 定 し た 各 流 星 の 実 経路 で , 前 報 と 同様 の 形式 で , 出 
現 点 及 び 消 滅 点 の 経度 , 緯度 及び 地上 の 高 さ (H), 経路 の 長 さ (①, 100 km の 距離 か 
5 見 た 標準 等 級 Cmo), AOD, 赤 紀 最 後に 所 属す る 流星 群 を 示し た 。 


oO 
Abbr. Observer Observed place L(E) ?(N) 
O x Jil BE (Onagawa Gr.) 宮城 県 女川 141°28/ 38°26/ 
永山 | Nagayama (Nm), 倉持 Kuramochi (Cr), 近藤 Kondo (Co) 其他 
S」 仙台 1 班 (Sendai Gr. I) 宮城 県 仙台 市 140 52 38 15 
高橋 仁 Z. Takahashi (Tk) 其他 
Ss 仙台 2 班 (Sendai Gr. IL) 同 a9 14052 38 15 
43325 Y. Konno (Ky), #fA4 M. Takeuchi (Tc), 助 野 Suxeno (Se) 
Mt x 戸 % (A. Mito) ‘Ee SUR EEA 14052 38 2 
Ab 阿部 H 3 CS. Abe) SST 138 32.4 37 21.6 
Mu # 山 @# = CY. Murayama) 同 = FAB 138 8.7 36.59.2 
Ty ms Wi A (CH. Takayama) pas BN 139 56.7 26 35.7 
Oy 3 野 f& # CY. Ogino) ERB AT 13919 3611 
KA 春日 部 高校 生 (CKasukabe Gr.) 
Cw WwW 井 ft (T. Kawai) 
Mj 宮地 義朗 CY. Miyaii) 同 春日 部 市 139 44.6 35 58.6 
Mo 桃 木 4e (H. Momoki) 
: 同 上 139 44.6 35 58.6 
本 Et Senzaitt) } A 谷塚 町 139 49 35 49 
No 長岡 信 i. CN. Nagaoka) 東京 都 北 区 上 十条 139 45 35 44 
M 武蔵 高校 生 (MusashiGr.) Al 練馬 区 139 40 35 44 
伴 場 次 郎 J. Bamba (Bb), 松 薬 谷 明 T. Matsuba (Mb), 吉川 省吾 S.Yoshikawa (Ys) 其他 
Iz fh HR 進 CS. Tizuka) 同 杉並 区 荻窪 18937.6 35 42.4 
Cg 熊谷 3H — (M. Kumagai) 同 濾 谷 区 代 * 木 139 40.7 35 40.1 
Kn & 36 Be #1 ME C1. Konko) 同 Gil 139 40.6 35 40.0 
N 日 本 学園 生 (Nihon Gakuen Gr.) 同 世 田 ヶ 谷 区 松原 町 139 39.1 35 39.8 
AVL Murayama (Ma) 其他 
Sk 坂巻 3 WE (S. Sakamaki) 同 = A RB 139 41.1 35 29.3 
Toe HH AW 5a (K. Takeuchi) 同 | (RARE 139 43.4 35 38.8 
Kr /) 32 i — (K. Kodaira) 神奈 川 県 川崎 市 139 37.5 35 33.6 
Ho #& SF 実 (M. Hoshino) 東京 都 町 田町 139 27 35 32 
So 佐藤 #b 夫 (K. Sato) AYR) LR CT 139 19.0. 35 23.5 
Co +R A= KR CF. Komatsu) 同 松田 町 139 18.5 35 20.5 
.K 鎌倉, 横浜 国 大 生 (Kamakura Gr.) lal 鎌倉 市 139 33.6 25 19.3 
“WPkGMI H. Kobayashi (Cb. 1953), 小林 謙一 郎 K. Kobayashi (Cy), 中 根 三 郎 S. Nakane 
(Nk) 其他 
Sn 庄野 ¥ # (CY. Shono) 同 大 磯 町 13918 35 18 


Ub- & % & FH (C. Ubukata) A | BP 139 34.9 35 17.4 
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In FF - # CY. Inouye) AHA R/S 139° 9/.7 35° 14/.8 
Il JAA, BORER Cwamura Gr.) lal ay 139° 8:7 35 9.4 
小林 部 叫 H. Kobayashi,(Cb. 1954), 小 出 康雄 Y. Koide (Cd), 近田 久美 子 K. Chikada 
(Ck) 其他 
Kd 神 A pe (S. Kanda) MAS SR 139 4.7 35 8.4 
Ns BE ¢ (CH. Nagafusa) 同 焼 準 市 138 19.4 34 51.7 
Cu ££ 林 良治 (Y.Kurebayashi) 同 金谷 町 138 8.7 34 48.6 
aa, 2. GW 
1950 s7T Obs Wt. Dur. Mag, Vel. Col. Begining Ending 
Aug. x 36 ee 
ioe aan) s m 
ofKe Sr 0.5 4 rR te peptee * ano 7 
ge 21-19 We NS 2 rS W 270 +29.5 253 +30.5 
Kn 5 0.5 3 M W 16 +34 22 +28 
9 1 43.9 on 4 0.6 3 Ss Ww 29 +48 46 +42 
‘Ker oS 0.3 4 M ia, 1654-62.5 $50 +68 
2 8.5 1g, Baa One 3 R BW 48 +60 62 +76 
TA ee = A = B 8+41 358.5433 
NST 5 0.1 1 vR BW 38 +56 33 +57 
(Kn 5 1.3 2.5 rs WwW 243.5—4 240 —21 
20 27 4G, a7 Se oe veS . YW 2660-11 955 235 
Kn ey a 3 Ss ~~ 2394-11 BE 1b39 49.5 
20 31.7 icp hy 0:5 3 rS rS 255 +37 240 +40? 
K 3 0.2 4 R — 258 + 8 258 — 1 
20-379 bi 0.1 1 R YY. 272-426 ** 275 +23 
ic 1 0.8 1.5 M Y? 258.5-—4 254 —11 
Aizen ¢ : 4 M ae OSes $250 38 
S 5. 0.6 0 Ss Ya 265410 967 {78 
1 a 2.0 0 S x 296 +0 280 —25 
91K A 2G 1 S 29205 ~975°— 380 
Be tee! yf 9 ee ee vS WY 809 +14 290,5—22 
5 0.7 5 M = 257 + 8 252 — 2 
21 17.1 {Se len 003 1 R BW 269.5426 161 +3 
123K x | =a 5 z69 —30 
ape {i 3 る =o S Ys - 312-24 268 —25 
0.4 4 M <4 2020540 COS 
oe fam 1!) tO Bw A om 2 
0.2 A M Ee 205045 oy 244 = 17 
BP OE Oe 3i- sR 0 BW 260'4¢10” “257 +.6 
0.3 2 R W 50 +46 37 +33 
14 14 a Od 3 R YR 61451 64 +43 
3 ? R WB 33 +62 26 +71 
147.5 1 nfl es 4 ey gl 1G Sie MN 
Iz 2 0.8 3 M W 28.5435 295 +26 
oa. 3 02 2 R W  23.5441.5 28 +29 
Se, of 04 4 M YR 49 +49 54 +33 
1.5 1-0 vR Wo =o 8 7 2 5 
2 89 {Kn ee ike Mb oh ah RT ES 
So oun 0.8 = vR WwW 38 +27 36 +21 
0.4 ein W 32.5421 27 +14 
ee ne Wa oh 4 R BW 50.5429 52.5419 
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2, & 2) 
: gini Endi 
No, Me «S.T. Obs. We. Dur. Mag Vel, Col, -Besining nding Rom 
Bem Ss m 
(Iz 4 0.7 3 rR W? = 56°-+.23° 56° +12° 
5019 14 2°20 1 co 3 03 2 rR BW 61 +28 65 +22 
No 5 0.5 2 R WwW 49 +19 51.5+32 T1.0s 
20 2 22 { So 40° 0.3 3 R B 59 +25 63.5+34 
Co 2 0.2 3 R YW C623 68 +32.5 
No 4 0.8 1 M WwW 30 +21 pe CY 
21 251 Cm ~ 3 1.2 2 rR Y 30 +21 22.5+12.5 T 2s 
So 5 0.4 1 R WB 47 +39 49.54+26.5 T 
No 5 0.9 06° OR W 8 +40 322 4374 CE 
22 2 58. of Kn 4 -- 1 M B 23 +40 353 +43 Ts 
So 5 0.5 1 vR WB 26 +60 5 6° 
23 es | 3 0.7 3 M Y 55 +44 57 +40 
So 3 0.4 251k WB 83 +46 89.5+40 
( Ho 4 0.8 1 rR W 55.5427 59 +22 T 
aA SSE CA CAT cm dang 0-8 1 ba — 7% +37.5 78431 T 
Kn 3 0.3 5 M — ~ 9725-10, 579516 
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1950 年 8 月 の も の は 日 本 天文 研究 会 に 報告 され た 流 性 観測 が , 観測 者 数 25 名 以上 , 
観測 回 数 113 夜 , 観測 時 数 延 176 時 間 , 観測 流星 数 1750 個 以 上 ある が , 未だ 整理 
され て いな か つた の で , それ を 調査 し て 32 個 の 実 経路 を 決定 する こと が で きた 。 実 経 
路 を 決定 し 得 た も の は すべ て 関東 地方 の 観測 者 の 観測 に よる も の で る ある る 。 

1952 年 8 月 の も の は , 前 報 の 原稿 を 完成 し た 後に 得 た 資料 に よ つ て 実 経路 を 決定 
得 た も の 13 Mle AACSB Lio ZO Nos. 5225~5229 の 5 (tae each ee 
会 員 の 観測 か ら 決 定 し た も の で あり , 観測 地 は 宮城 県 女川 , MA, HED 3 個所 で あ 
Do COMO 8 個 は 静岡 県 , 長野 県 並 に 関東 地方 の 上 冠 の も の で ある 。 

1953 年 8 月 の も の は , 前 報 の も の と ほ 革 同じ く 関 東 地方 , 静岡 県 , 新潟 県 の 観測 者 
の 観測 に よる も の で , 横浜 国立 大 学 学芸 学部 天 六 気 祭 部 学生 並 に 日 本 天文 研究 会 会 員 
の 観測 で ある 。 

1953 年 8 月 11 日 Oh 138m~3h 11m N.S. T. の 間 , 埼玉 県 春日 部 高校 天象 部 の 流星 
観測 に 参加 し た 同 校 卒業 生 妊 崎 秀 夫 I. Susaki 氏 が コ ニ = カエ へ キサ ー『F 3.5 FEB 
に より 撮影 し た 写真 を カビ ネ 大 に 引 仲 し た 印画 を 途 ら れ た 。 HERO MHRA RO 
ダグ ルー プ ・ カ ウト の 記録 に よれ ば 3h 6m 56s, 東京 の 坂 大 氏 に よれ ば 3h7m.4 と 記録 さ 
れ て いる も の で を る り , 印画 の 上 か ら 筆 者 の 決定 し た 発光 点 a 48.0, 54+23°15/, 消 減 
点 4h17m.9, 84-14°32! (1953.0) な る 値 を 得 た 。 こ の 概 炎 は 第 2 HO No. 5323 の 
中 KA と 示し た 項 に 記し た 。 WM UAOE, SRR, 色 は ダル ー プ ・ カ ウッ ント の 思 
a eR 

1954 年 8 月 の 7 個 は すべ て 横浜 国立 大 学 学芸 学部 天文 気 祭 部 学生 の 観測 に よる る も 
OC, 鎌倉 の 学芸 学部 と , 新設 の 神奈 川 県 足柄 下 郡 岩村 の 理科 教育 実験 所 と で 行わ れ 
RE 

実 現 観測 か ら 求 め た 実 経路 で ある か ら , 第 3 実 の 数 値 は か な り 不 正確 の も の が あろ 

と 思わ れる 。 消滅 点 の 高 さ H CMR +t の 次 に 示し た 数 値 は 平均 偏差 で これ は 経 
路 の 精 鹿 を 示す 資料 と し て 示し た 。 + DOU て な いも の は 消滅 点 の 高 さ と し て 一 施 の 
観測 者 の 結果 だ け を 探 用 し た 場合 で ある 。 RHR RR.P. に 下線 を つけ た も の は , B 
測 不良 の た め , その 数 値 を 仮定 し た 事 を 示す 。 

前 報 の 結果 と 合せ て 103 個 は すべ て 8 月 8~15 日 の 観測 で あり , “reo ame 
56 A, ヵ カシ ォ ペ ィ イプ アプ 座 群 13 個 , MERE 7 個 で , その 他 は 3 個 以 下 で ある 。 HAR 
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が 増し た 結果 で あろ うと 選 わ れる 。 。 6a 
終り に 多数 の 観測 を 報告 され た 方 々 と , Afe MD OH DBI ERR 
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